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Abstract. We use milli-arcsecond accuracy astrometry 
(proper motions and parallaxes) from Hipparcos and from 
radio observations to retrace the orbits of 56 runaway stars 
and nine compact objects with distances less than 700 
pc, to identify the parent stellar group. It is possible to 
deduce the specific formation scenario with near certainty 
for two cases, (i) We find that the runaway star C Ophiuchi 
and the pulsar PSR J1932+1059 originated about 1 Myr 
ago in a supernova explosion in a binary in the Upper 
Scorpius subgroup of the Sco 0B2 association. The pulsar 
received a kick velocity of ^350 km s~^ in this event, which 
dissociated the binary, and gave C Oph its large space 
velocity, (ii) Blaauw & Morgan and Gies & Bolton already 
postulated a common origin for the runaway-pair AE Aur 
and fj. Col, possibly involving the massive highly-eccentric 
binary i Ori, based on their equal and opposite velocities. 
We demonstrate that these three objects indeed occupied 
a very small volume ~2.5 Myr ago, and show that they 
were ejected from the nascent Trapezium cluster. 

We identify the parent group for two more pulsars: 
both likely originate in the ^ 50 Myr old association Per 
0B3, which contains the open cluster a Persei. At least 
21 of the 56 runaway stars in our sample can be linked to 
the nearby associations and young open clusters. These in- 
clude the classical runaways 53 Arietis (Ori OBI), ^ Persei 
(Per 0B2), and A Cephei (Cep 0B3), and fifteen new iden- 
tifications, amongst which a pair of stars running away in 
opposite directions from the region containing the A Ori 
cluster. Other currently nearby runaways and pulsars orig- 
inated beyond 700 pc, where our knowledge of the parent 
groups is very incomplete. 

Key words: Astrometry - Stars: early-type - Stars: in- 
dividual: Runaway ~ Stars: kinematics - Pulsars: general 
- Supernova: general 



1. Introduction 

About 10-30% of the O stars and 5-10% of the B stars 



(Gies 1987; Stone 1991) have large peculiar velocities (up 
to 200 km s~^), and are often found in isolated locations; 
these are the so-called 'runaway stars' (Blaauw 1961, here- 
after Paper I). The velocity dispersion of the population 
of runaway stars, cr^ ^ 30 km s~^ (e.g., Stone 1991), is 
much larger than that of the 'normal' early-type stars, 
ay ~ 10 km s""'^. Besides their peculiar kinematics, run- 
away stars are also distinguished from the normal early- 
type stars by an almost complete absence of multiplicity 
(cf. the binary fraction of normal early-type stars is >50% 
[e.g.. Mason et al. 1998| ). Furthermore, over 50% of the 



(massive) runaways have large rotational velocities and 
enhanced surface helium abundances (Blaauw 1993). 

Several mechanisms have been suggested for the origin 
of run away stars (Zwicky |l957 ; Paper I; Poveda, Ruiz fc 
Allen |l967t Carrasco et al. |l980t Gies & Bohon |1986D , 
two of which are still viable: the binary-supernova scenario 
(Paper I) and the dynamical ejection scenario (Poveda et 
al. 1967). We summarize them in turn. 



Binary- supernova scenario (BSS): In this scenario a run- 
away star receives its velocity when the primary compo- 
nent of a massive binary system explodes as a supernova. 
When the supernova shell passes the secondary the grav- 
itational attraction of the primary reduces considerably, 
and the secondary starts to move through space with a 
velocity comparable to its original orbital velocity (30- 
150 km s^^). What remains of the primary after the explo- 
sion is a compact object, either a neutron star or a black 
hole. Depending on the details of the preceding binary 
evolution, the eccentricity of the orbit, and the kick veloc- 
ity fkick due to the asymmetry of the supernova explosion 



(e.g.. Burrows, Hayes & Fryxell 1995), the compact object 
may or may not remain bound to the runaway star (e.g., 



Hills 1983). In most cases less than half of the total mass 
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of the binary is expelled in the explosion. As this is insuffi- 
cient to break up the binary (Paper I), most BSS runaways 
are expected to remain double. Examples of such systems 
are provided by the high- mass-X-ray binaries (e.g., van 
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den Heuvel et al. 2000). Their typical velocities of ^50 
km (Kaper et al. |1997t Chevalier & Ilovaisky [1998D 
are the natural result of the recoil velocity acquired when 
the supernova shell leaves the binary system. 

Several searches failed to find compact companions of 
classical runaway stars (Gies & Bolton 1986; Philp et al. 
1996 ; Sayer, Nice & Kaspi 1996| ), suggesting that for these 



systems I'kick of the neutron star was large enough (several 
100 km s~^) to unbind the binary (e.g., Frail & Kulkarni 
1991|; Cordes, Romani & Lundgren |1993| ; Lai [19991 ). The 



certain (e.g., Hills 1983; Lorimer, Bailes & Harrison 1997; 
Hansen & Phinney p"997[ Hartman |1997D . 

Single BSS runaways must originate in close binaries 
because these systems have the largest orbital velocities, 
and therefore they have experienced close binary evolu- 
tion before being ejected as a runaway. This leads to the 
following observable characteristics: 

BSS runaways are expected to have increased helium 
abundance]^ and large rotational velocity: when the pri- 
mary fills its Roche lobe, mass and angular momen- 
tum is transferred to the star that will become the 
runaway. The mass transfer stops when the 'primary' 
has become a helium star (i.e., only the helium core 
remains). This process enriches the runaway with he- 
lium, and spins it up (e.g.. Packet 1981; van den Heuvel 



1: 



1985; Blaauw 1993) 



2: A BSS runaway can become a blue straggler, because it 
is rejuvenated during the mass-transfer period through 
the fresh fuel it receives from the primary. 

3: The kinematic age of a BSS runaway star (defined as 
the time since the runaway left its parent group) should 
be smaller than the age of the parent group. The pri- 
mary of the original binary system first evolves for sev- 
eral Myr before it explodes and the runaway is ejected. 

Dynamical ejection scenario (DES): In this scenario run- 
away stars are formed through gravitational interac- 
tions between stars in dense, compact clusters. Although 
binary-single star encounters produce runaways (e.g.. Hut 
& Bahcall 1983), the most efficient intera ction is the en- 
counter of two hard binary systems (Hoffer 198S ) . Detailed 
simulations show that these collisions pro duce runaways 
with velocities up t o 200 km s~ ^ (Mikko la |l9834 |l983b ; 
Leonard & Duncan |l984 |l990| ; Leonard |l991^ . The out- 
come of a binary-binary collision can be (i) two binaries, 
(ii) one single star and a hierarchical triple system, (iii) 
two singl e stars and one binary, and (iv) four single stars 
(Leonard 1989 ). In most cases the collision will result in 
the ejection of two single stars and one hard binary with 
an eccentric orbit (Hoffer 1983; Mikkola 1983a ). Since the 
resulting binary is the most massive end product of the 



^ The abundances of other elements are increased as well, in 
particular that of nitrogen which can be enhanced by a factor 
4-5 (Van Rensbergen priv. comm.) 



collision it is unlikely to gain a lot of speed; it might even 
remain within the parent cluster. This process naturally 
leads to a low (0-33%) runaway binary fraction which is in 
qualitative agreement with the observations (e.g, Gies & 
Bolton 1986). For the DES to be efficient the initial binary 
fraction in clusters needs to be large. Recent observations 
show that the binary fraction for massive stars in young 
clust ers is >50% (Abt |1983|; K roupa, Petr & McCaugh- 
rean |1999| ; Preibisch et al. [1999| [approaching 100%]). 



typical magnitude of this 'threshold' kick velocity is un- \- 



DES runaways have the following characteristics: 

DES runaways are formed most efficiently in a high- 
density environment, e.g., in young open clusters. They 
may also originate in OB associations. These birth 
sites of massive stars are unbound stellar groups and 
therefore expand (e.g., Blaauw 1952a, 1978 Elmegreen 



1983 ; Kroupa 2000a), so DES runaways must have been 



ejected very soon after their formation. The kinematic 
age and the age of the parent association are thus 
nearly equal. 

DES runaways are not expected to show signs of bi- 
nary evolution such as large rotational velocities and 



increased helium abundance. However, Leonard (1995) 



suggested that some binary-binary encounters produce 
runaways consisting of two stars that merged during 
the interaction. These would have enhanced helium 
abundances and large rotational velocities (Benz & 
Hills 1987, but see Lombardi, Rasio & Shapiro 1995). 



3; DES runaways are expected to be mostly single stars. 

Which of the two formation processes is responsible for 
runaway stars has been debated vigorously. Both mecha- 
nisms create stars with large peculiar velocities which en- 
able them to travel far from their parent group: a velocity 
of 100 km s~^ corresponds to '--^lOO pc in only 1 Myr. 
The relative importance of the two scenarios can be es- 
tablished by (i) studying the statistical properties of the 
ensemble of runaway stars, or by (ii) investigating indi- 
vidual runaways in detail. The former approach is based 
on differences in the general runaway characteristics pre- 
dicted by each scenario. This requires a large, complete 
database of runaway stars, which is, to date, unavailable 
(e.g., Moffat et al. 1998). Here we therefore follow the lat- 
ter, individual approach by retracing the orbits of runaway 
stars back in time. The objects encountered by a runaway 
along its path (e.g., an open cluster, an association, other 
runaways, or a neutron star), and the times at which these 
encounters occurred, provide information about its forma- 
tion. Evidence for the BSS as the formation mechanism for 
single runaway stars is to find a runaway and a neutron 
star (pulsar) which occupied the same region of space at 
the same time in the past. Evidence for the DES is to find 
a common site of origin for the individual components of 
the encounter, e.g., a pair of runaways and a binary, in a 
dense star cluster. 

The individual approach requires highly accurate posi- 
tions (a,(5,7r) and velocities (/Xc(*,M5,Vrad)- Here a denotes 
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Table 1. The nearby runaway stars and pulsars with accurate astrometry. HIP indicates the number of the runaway 
star in the Hipparcos Catalogue, Wspacc indicates the space motion of the runaway star relative to Galactic rotation 
(in km s~^), and PSR indicates the pulsar identifications. 



HIP 




HIP 


"f^spacc 


HIP 


t'space 


HIP 


f space 


HIP 


'^^spacc 


HIP 


"l^spacc 


PSR 


3478 


80.4 


28756 


196.6 


43158 


57.2 


61602 


30.2 


91599 


44.7 


101350 


36.4 


J0826+2637 


3881 


32.1 


29678 


63.0 


45563 


125.9 


62322 


43.9 


92609 


31.0 


102274 


46.1 


J0835-4510 


9549 


107.9 


30143 


55.5 


46928 


45.3 


66524 


112.7 


94899 


162.8 


103206 


32.3 


J0953+0755 


10849 


50.0 


35951 


34.9 


46950 


32.1 


69491 


77.2 


94934 


94.0 


105811 


38.3 


J1115+5030 


14514 


39.4 


36246 


32.1 


48715 


34.5 


70574 


205.3 


95818 


34.7 


106620 


45.2 


J1136+1551 


18614 


64.9 


38455 


41.4 


48943 


35.2 


76013 


69.0 


96115 


165.6 


109556 


74.0 


J1239+2453 


20330 


34.7 


38518 


31.1 


49934 


31.2 


81377 


23.5 


97774 


35.0 






J1456-6843 


22061 


86.5 


39429 


62.4 


52161 


34.8 


82171 


62.9 


97845 


70.3 






J1932+1059 


24575 


113.3 


40341 


61.7 


57669 


31.1 


82868 


30.3 


99435 


39.4 






Geminga 


27204 


107.8 


42038 


31.3 


59607 


78.8 


86768 


30.1 


99580 


55.6 









right ascension, S declination, tt parallax, fia* — HaCosS 
proper motion in right ascension, fis proper motion in dec- 
lination, and Wrad the radial velocity. The milli-arcsecond 
(mas) accuracy of Hipparcos astrometry (ESA 1997) al- 
lows specific investigations of the runaway stars within 
'--^TOO pc. Positions and proper motions of similar accu- 
racy are now available as well for some pulsars through 
timing measurements a nd VL BI observa tions (e.g., Tay- 
lor, Manchester & Lyne |l993| ; Campbell |l995|) . The Hip- 
parcos data also significantly improved and extended the 
mem bershi p lists of the nearby OB associations (de Zeeuw 
et al. 1999|) , and of some nearby young open clusters. The 
resulting improved distances and space velocities of these 
stellar aggregates make it possible to connect the run- 
aways and pulsars to their parent group, and, in some 
cases, to identify the specific formation scenario (Hooger- 
werf, de Bruijne fc de Zeeuw 2000| ; de Zeeuw, Hoogerw- 
erf & de B ruijne ^000 ). Pre-Hipparcos data (e.g., Blaauw 
& Morgan 1954 ; Paper I; Bl aauw 1993| ; van Rensbergen, 
Vanbeveren & de Loore 1996) allowed identification of the 
parent groups for some runaways, but generally lacked the 
accuracy to study the orbits of the runaways in detail (but 
see Blaauw & Morgan |1954| and Gies & Bolton |1986| ). 

We define a sample of nearby runaways and pulsars 
with good astrometry in and then analyse two cases 
in depth: C Oph and PSR J1932-hl059 in §|and AE Aur, 
H Col and l Ori in We apply the method developed 
in these sections to the entire sample of runaways and 
pulsars in and §0. We discuss helium abundances, 

rotational velocities and the blue straggler nature of run- 
aways in §H and §^ and summarize our conclusions in § |lO| . 

2. Nearby runaway stars and pulsars 

The parent group is known for about a dozen 'classical' 
runaway stars (Paper I; Blaauw 1993). The Hipparcos 



dates in Hipparcos Proposal 141q (de Zeeuw et al. 1999). 
Many of these objects are located beyond ~700 pc, where 
the Hipparcos parallax measurement is of modest qual- 
ity. For this reason we restricted ourselves to a sample of 
nearby runaway stars, and added to this the (few) nearby 
pulsars with measured proper motions. 

2.1. Selection of the sample 

We started with all 1118 O to B5 stars in the Hipparcos 
Catalogue which have radial velociti es list ed in the Hip- 
parcos Input Catalogue (Turon et al. 1992 ). Next we only 
considered those stars which have significant parallaxes 
(tt — 2cr^ > mas) and proper motions (ct^/a* < 0.1), and 
space velocities larger than 30 km s~^ with respect to the 
standard of rest of the runaway. For the last requirement 
we corrected the runaway velocity for Solar motion and 



Galactic rotation (Dehnen & Binney 1997). The somewhat 



arbitrary choice of the velocity limit of 30 km s~ mini- 
mizes the contamination of the sample by normal O and B 
stars (§0). These criteria yield 54 runaway candidates (five 
of which are classical runaways)]^ This new sample does 
not contain the nearby runaways C Oph and ^ Per. The 
former is not selected because its space velocity is smaller 
than 30 km s^^ (although its velocity relative to its par- 
ent group Sco 0B2 is larger, cf. §||) and the latter is not 
selected because cr^//x > 0.1. However, since the runaway 
nature of these two stars is well established (e.g., Paper 
I) we included them in our sample, bringing the total to 
56. The Hipparcos numbers and space velocities of these 
56 stars are listed in Table ^ Panel a of Figure |l| shows 
the histogram of the derived space velocities. 



Catalogue contains these stars, as well as many additional 
O and B stars which were known in 1982 to have large 
radial velocities, including 153 of the 162 runaway candi- 



^ From a list of O and B stars with Wrad > 30 km s^^ provided 
by the late Jan van Paradijs. 

^ The sample criteria excl ude the known runaways, e.g., 
a Cam, 68 Cyg (Blaauw [l993| ), HIP 35707 (Co mero n, Torra & 



Gomez |1998| ), and HIP 60350 (Maitzen et al. |l998[ ) which aU 
have distances larger than 1 kpc. 
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Fig. 1. Panel a: Histogram of the spac e motions of the sample of runaway stars defined in §2.1, Panel b: Distribution 
of pulsars from the Taylor et al. (1993) catalogue with measured proper motions. The light grey histogram shows all 
pulsars within 2 kpc and the dark grey histogram shows the pulsar with accurate proper motions {cr^j./ H < 0.1). The 
latter, for D < 1 kpc, is the pulsar sample defined in §2.1, 



We selected a sample of nearby pulsars from the Tay- 
lor, Manchester & Lyne ( 1993|) ca talogue, as updated on 
http : //pulsar . princeton . edu/l . It contains 94 pulsars 
with known proper motions and distances. Only seven of 
these meet our distance {D < 1 kpc) and proper motion 
{(Ji^/li < 0.1) constraints (see panel b of Figure |l|). Most 
pulsar distances are derived from the dispersion measure. 
These distances are unreliable, especially for nearby ob- 
jects, since they depend on the local properties of the ISM. 
For one nearby pulsar, PSR J0953-f0755, high precision 
VLBA mea surements became available recently (Brisken 
et al. 2000 ). We added this pulsar to our sample. The 
eight pulsars are listed in Table |l|, together with Geminga, 
a nearby neutron star which is not a pulsar, for which an 
accurate proper motion is known (Caraveo et al. 



1996) 



Our sample of nearby runaway stars and compact ob- 
jects is severely incomplete. The Hipparcos Catalogue is 
complete to V = 7.3-9 mag, with the limit depending 
on Galactic latitude and spectral type (2163 of the 3622 
O to B5 stars have V > 7.3 mag). The data available 
for the O and B stars is inhomogeneous and incomplete, 
e.g., less than a third of the O to B5 stars in the Catalog 
has a measured radial velocity. We have excluded those 
with large v^ad but insignificant proper motions, as their 
retraced orbits are uncertain. The beamed nature of the 
radio emission from pulsars hides many from observation, 
and not all of those that do radiate in our direction have 
been found. Of these, only a few have an accurately mea- 
sured proper motion and a reliable distance. 

2.2. Nearby OB associations and open clusters 

We adopt the positions and mean space motions of the 
OB associations with in 700 pc of the Sun as derived by 
de Zeeuw et al. ( 1999 ) from Hipparcos measurements. For 



the open clusters we compiled a list from the WEBDA cat- 
alogue ( http : / / obswww . unige . ch/webda/| ) , and consider 
only those which are young (r < 50 Myr) and with dis- 
tances less than 700 pc as likely parent groups. The age 
requirement is comparable to the age of the oldest run- 
aways we consider here (B5V). Typical pulsar ages are 
less than 50 Myr (e.g., Blaauw & Ramachandran 1998) 



This selection yields nineteen open clusters (see Table ^), 
of which five are already covered in the study of the nearby 
associations by de Zeeuw et al. (1999). To obtain the space 
motion of these clusters we use the WEBDA member stars 
listed in the Hipparcos Catalogue to obtain reliable as- 
trometry, and those in the Hipparcos Input Catalogue to 
obtain the radial velocity. In this way we are able to con- 
struct a more or less reliable space motion for seven of 
the fourteen remaining open clusters (those labeled "Y" 
or "?" in Table | which summarizes the results). 

2.3. Orbits 

Traditionally, the orbits of runaway stars have been traced 
back in time using straight lines through space. This is 
sufficiently accurate for identification of the parent group 
for times up to a few Myr and distances less than a few 
hundred pc. To make sure we include the effect of the 
Galactic potential, we use a fourth-order Runge-Kutta 
numerical integration method, with a fixed time-step of 
10 000 yr, to calculate the orbit. The Galactic potential 
we use consists of (i) a logarithmic potential for the halo, 
(ii) a Miyamoto-Nagai potential for the disk, and (iii) a 
Plummer potential for the bulge of the Galaxy. The po- 
tential predicts Oort constants A = 13.5 km s"^ kpc~^ 
and B = —12.4 km s^^ kpc^^ and a circular velocity 
i^circ — 219.8 km s~^ at Rq = 8.5 kpc. These values agree 



with those which Feast & Whitelock ( 1997 ) obtained us 
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Table 2. The nearby {D < 700 pc), young (r < 50 Myr), open clusters listed in the WEBDA Catalogue. The table gives 
the designation of the open cluster (Name), its position on the sky (£, b), the number of member stars contained in the 
Hipparcos (HIP) and Hipparcos Input Catalogues (HIC), its distance, proper motion, and radial velocity as obtained 
from the Hipparcos data {D, [/i^,, ^j,], Wrad), and whether the cluster is a candidate parent group in this study (Cand.). 
The candidate status is denoted by "Y" for the clusters with well-determined positions and velocities, by "N" for 
clusters for which not all information is available (either astrometry or radial velocity), by "N*" if the astrometry does 
not show a clear signature of an open cluster (i.e., a clump in the proper-motion vs. proper-motion diagram), or by 
"?" if the measurements are not very reliable (either because of a small number of member stars or a large spread in 
the data). 



Name 


I 

[deg.] 


b 

[deg.] 


HIP 

[#] 


HIC D fl£, fib t'rad 

[#] [pc] [mas yr^-""] [mas yr^"""] [km s^"""] 


Cand. 


CoUinder 359 


29.75 


12.54 


10 




N* 


IC 4665 


30.61 


17.08 


13 


5 385 ±40 -7.2 ±0.3 -3.0 ± 0.3 -13.5 ±3.0 


Y 


Stephenson 1 


66.85 


15.51 







N 


Roslund 5 


71.40 


0.25 


13 




N* 


Stock 7 


134.68 


0.04 


3 




N* 


a Persei 




Central part of the Per OB3 association, contained in the list of de Zeeuw et al. 




IC 0348 








Associated with Per OB2 




CoUinder 69 


195.05 


-12.00 


6 


— A Ori cluster 


N* 


NGC 1976 








Trapezium cluster: associated with Ori OBI 




NGC 2232 


214.36 


-47.65 


10 


3 365 ±40 0.7 ±0.5 -5.2 ± 0.5 14.6 ± 3.0 


? 


CoUinder 121 






Contained in list of nearby associations of de Zeeuw et al. 




CoUinder 140 


245.18 


-7.87 


14 


4 375 ±40 -7.4 ±0.5 -5.5 ± 0.5 22.4 ± 3.0 


Y 


CoUinder 135 


248.76 


-11.20 


19 


4 300 ±30 -10.3 ±0.5 -6.8 ± 0.5 16.4 ± 3.0 


Y 


Pismis 5 


259.39 


0.86 







N 


Pismis 4 


262.74 


-2.37 


4 





N 


Trumpler 10 






Contained in list of nearby associations of de Zeeuw et al. 




IC 2391 


270.36 


-6.88 


24 


13 150 ± 30 -33.1 ±0.5 -6.0 ± 0.5 15.0 ± 3.0 


Y 


vdB-Hagen 99 


286.56 


-0.63 


7 


2 500 ±50 -13.1 ±0.5 -6.4 ±0.5 12.0 ± 3.0 


? 


IC 2602 


289.60 


-4.90 


25 


8 140 ±10 -20.4 ±0.5 1.2 ± 0.5 24.1 ±3.0 


Y 



ing Hipparcos data: A ~ 14.82 ±0.84 km s ^ kpc ^, B — 2.4- Identification of parent groups 
-12.37±0.64kms-i kpc-\ Wcirc = 231.2±16.2 km s'^ at 
Ro = 8.5 kpc. Since the volume covered in the orbit inte- 
gration is typically a few hundred pc, and the time of the 
integration is typically less than 10 Myr, perturbations of 
the orbits caused by small-scale structure in the disk are 
negligible. 

Before integrating the orbit, we correct the observed 
velocity f, for (i) the Solar motion with respect to the 
Local Standard of Rest, v\sr (Dehnen & Binney 1997), and 
(ii) the Galactic rotational velocity of the Local Standard 



of Rest, (Binney & Tremaine 1987, p. 14). The stellar 



velocity Vgai relative to the Galactic reference frame is 
then given by 



(1) 



To retrace the orbit, we reverse the velocity and integrate 
forward in time. We calculate the distance of a star as 
I/tt, where tt is the trigonometric parallax. Since we use 
the individual parallax, we cannot correct this distance for 



possible biases (e.g.. Smith & Eichhorn 1996) 



We calculate the past orbit of each of the 56 runaway stars 
listed in Table |l] for 10 Myr. We do this 10 000 times for 
each star, in order to sample the error ellipsoid of the mea- 
sured parameters, defined by the covariance matrix of the 
Hipparcos astrometry and the error in the radial velocity 
measurement. Retracing the orbit of a pulsar is more diffi- 
cult, because the radial velocity is unknown. We therefore 
cover a range of radial velocities of Wiad = ± 500 km 
in the orbit integrations for the pulsars. Figure ^ shows the 
positions of the runaways and pulsars on the sky, together 
with their orbits, retraced back for only 2 Myr so as not to 
confuse the diagram. Three orbits are shown for each pul- 
sar: for Vrad = km s^^ (filled square), Vrad = 200 km s^^ 
(open square), and Wrad = —200 km s^^ (open star). 

We also retrace the orbits of the set of nearby OB as- 
sociations and open clusters defined in These groups 
have typical linear dimensions of 10-30 pc. We consider a 
group to be a possible site of origin for a runaway or pulsar 
if the minimum separation between the runaway/pulsar 
and the group was less than 10 pc at some time in the past 
10 Myr. With this definition, we find a parent group for 21 
of the 56 runaways. These stars are indicated by the filled 



6 



R. Hoogerwerf et al. 
Runaways 



On the origin of the O and B-type stars with high velocities II 
^- = 90° 




L = -90° 
I- = 90° 




L- = 90° 



Fig. 2. Top: Sample of runaway stars defined in |2Jj, in Galactic cordinates. The open circles denote the present 
positions of the runaways, and the arcs show their past orbits, calculated for 2 Myr. The filled circles are the runaways 
for which we can identify the parent association. The numbers refer to the entries in Table |3|. The asterisks indicate 
two additional runaways (72 Col, HIP 94899 [left most of the two asterisks]) discussed in §U The grey fields outline 
the nearby OB associations (de Zeeuw et al. 1999). From left to right and from top to bottom: Per 0B3 (a Persei), 
Per 0B2, Cep 0B3, Cep 0B2, Cep 0B6, Lac OBI, Upper Scorpius, Upper Centaurus Lupus, Lower Centaurus Crux, 
Tr 10, Vel 0B2, Col 121, and Ori OBI. The open clusters are identified by the filled light-grey circles for those 
with reliable positions and velocities, and by the open, crossed circles for the remaining clusters. The positions and 
designations of the clusters can be found in Table ^. Bottom: Pulsar sample defined in |2.l| , in Galactic coordinates. 
The filled circles indicate the present positions of the pulsars. The past orbits of pulsars, calculated for 2 Myr, are 
shown for three different assumed radial velocities: km s~^ (filled squares), 200 km (open squares), —200 km s~^ 
(open stars). The pulsars are labeled 1 through 8; 1: J0826+2637, 2: J0835-4510 (Vela pulsar), 3: J0953+0755, 4: 
J1115+5030, 5: J1136+1551, 6: J1239+2453, 7: J1456-6843, 8: J1932+1059. Number 9 is the neutron star Geminga. 
The associations and open clusters typically move comparatively little in 2 Myr. 

circles in Figure |, and include the seven classical run- PSR J0835-4510, PSR J1115-f5030, PSR J1239-I-2453, 
aways in the sample. We discuss them in detail in §§||-|| PSR J1932-I-1059, and Geminga (objects 1, 2, 4, 6, 8, and 
below. Six of the nine neutron stars possibly traversed one 9 in Figure |^) . We discuss them in ^ and §^1 and identify 
of the nearby stellar groups; these are PSR J0826-I-2637, the parent group for four of them. Table y summarizes 
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the data for the 22 runaways, four pulsars, and Geminga. 
The pulsars and runaways for which we cannot identify a 
parent group are discussed further in §|^. 



Table 3. Data for the nearby runaway stars and pulsars discussed in this paper. The stars which have an * appended to their HIP identifier are the classical 
runaways (i.e., a parent group was already known before this study, see Paper I. Unless indicated otherwise, the position (a, (5), proper motion (/ia*, /ia), 
and parallax (tt) were taken from the Hipparcos Catalogue (ESA 1997), the radial velocity (vrad) from the Hipparcos Input Catalogue (Turon et al. 1992), 
the space velocity (tigpace) with respect to the standard of rest of the runaway, the rotational velocity {vrot sini) for the runaways from Penny (1996) and 
the period P for the pulsars in seconds, the spectral type from Mason et al. (1998) or the Hipparcos Catalogue for the nmaways and the characteristic age 
(r = P/(2P)) for the pulsars, and the helium abundance (e) from Herrero et al. (1992), defined as the number of He atoms relative to H. The mass M^k 
has been derived from the Schmidt-Kaler (1982) calibration, using interpolation. The mass Mbb is taken from Vanbeveren, Van Rensbergen & de Loore 
(1998). The last column (N) indicates the number of the nmaway/pulsar in Figure 2. The proper motion and radial velocity are not corrected for Solar 
motion and Galactic rotation. The astrometric data (a, 5, tt, Ha*, and us) for the pulsars (the last five lines) are taken from the Taylor, Manchester & 
Lyne (1993) catalogue. Abbreviations used: mas = milli-arcsec; fia* = fJ-a cos S. 



HIP 


HD 


Name 


a (" 

[J1991.25] 


S (° ' ") 
[J1991,25] 


TT 

[mas] 


[mas yr~^] 


[mas yr~^] 


[km s~ 






[km s-i] 


Vrat sin i 

[km s-i; 


SpT 

1 


MsK 

[Mq] 


Mbb 


e 
[#] 


N 


3881 


4727 


L/ And 





49 


48. 


.83 


+41 


04 


44,2 


4,80 ± 0,75 


22.68 + 0.53 


-18.05 + 0.48 


-23.9 + 


1 


.2 


32.1 


80" 


B5V+F8V 


6.9" 






1 


14514* 


19374 


53 Ari 


3 


07 


25, 


.69 


+ 17 


52 


47.9 


4.32 ± 0.98 


-23.54 ± 0.93 


9.30 + 0.95 


21.2 + 


1 


.2" 


39.4 


10"* 


B1.5V 


10.4 


8.5 




2 


18614* 


24912 


5 Per 


3 


58 


57, 


.90 


+35 


47 


27.7 


1.84 ± 0.70 


1.92 ± 0.74 


2.30 + 0.62 


58.8 ± 


5 


.a" 


64.9 


204 


07.5III 


33.8 


33.5 


0.18 


3 


22061 


30112 




4 


44 


42 


.16 


+0 


34 


05.4 


2.94 ± 0.86 


-44.89 ± 0.77 


-29.28 ± 0.67 


6.0 ± 


5 


.0 


86.5 




B2.5V 


8.6 


7.5 




4 


24575* 


34078 


AE Aur 


5 


16 


18 


.15 


+34 


18 


44.0 


2.24 ± 0.74 


-4.05 ± 0.66 


43.22 + 0.44 


57.5 ± 


1 


.2 


113.3 


25 


09.5V 


15.9 


21.1 


0.09 


5 


26241 


37043 


L Ori 


5 


35 


25 


.98 


-05 


54 


35.6 


2.46 ± 0.77 


2.27 ± 0.65 


— 0.62 + 0.47 


28.7 ± 


1 


.if 


8.0 


71' 


09III+B1III'* 


37.8^ 


38.6 






27204* 


38666 


ti Col 


5 


45 


59 


.89 


-32 


18 


23.0 


2.52 ± 0.55 


3.01 ± 0.52 


-22.62 + 0.50 


109.0 ± 


2 


. .1 


107.8 


111 


09.5V 


15.9 


21.1 




6 


29678 


43112 




6 


15 


08, 


.46 


+ 13 


51 


03.9 


2.38 ± 0.72 


24.21 ± 0.76 


10.65 + 0.49 


36.0 ± 


5 


.0 


63.0 


<25^ 


BIV 


11.5 


12.0 




7 


38455 


64503 




7 


52 


38, 


.65 


-38 


51 


46.2 


5.09 ± 0,52 


-9.49 ± 0.43 


4.02 + 0.42 


-31.0 ± 


5, 


.0 


41.4 


212*= 


B2V 


9.4 


8.0 




8 


38518 


64760 




7 


53 


18, 


.16 


-48 


06 


10,6 


1.68 ± 0,50 


-4.90 ± 0.53 


5.89 + 0.38 


41.0 + 


5, 


.0 


31.1 


220-^ 


B0.5Iab 


25.0 


35.1 




9 


39429 


66811 


C Pup 


8 


03 


35, 


.07 


-40 


00 


11,5 


2,33 ± 0,51 


-30.82 ± 0.44 


16.77 + 0.41 


-23.9 + 


1 


.2 


62.4 


203 


04I 




67.5 


0.14' 


10 


42038 


73105 




8 


34 


09 


.60 


-53 


04 


17,5 


2,87 + 0,47 


-12.14 + 0.54 


10.13 + 0.48 


37.0 + 


10 


.0 


31.3 




B3V 


7.9 


7.0 




11 


46950 


83058 




9 


34 


08, 


.80 


-51 


15 


19,0 


3.50 ± 0,53 


-8.50 + 0.49 


6.39 + 0.48 


35.0 + 


10 


.0 


32.1 




B1.5IV 


10.4'" 


9.0 




12 


48943 


86612 




9 


59 


06 


.32 


-23 


57 


02,8 


5.19 ± 0.77 


-23.22 ± 0.70 


5.30 + 0.70 


39.0 ± 


5, 


.0 


35.2 


230"* 


B5V 


5.8 






13 


49934 


88661 




10 


11 


46, 


.47 


-58 


03 


38.0 


2.52 ± 0.50 


-10.71 ± 0.49 


6.63 + 0.45 


31.0 ± 


10 


.0 


31.2 


280"* 


B2IVnpc 


9.4'" 


8.0 




14 


57669 


102776 




11 


49 


41 


.09 


-63 


47 


18.6 


7.10 ± 0.69 


-17.93 ± 0.95 


4.44 + 0.63 


29.0 ± 


2 


.5 


31.1 


251" 


B3V 


7.9 


7.0 




15 


69491 


124195 




14 


13 


39, 


.84 


-54 


37 


32.2 


2.96 ± 0.63 


-18.03 ± 0.41 


-11.15 + 0.41 


66.0 ± 


10 


.0 


77.2 




B5V 


5.8 






16 


76013 


137387 


Aps 


15 


31 


30, 


.82 


-73 


23 


22.4 


3.20 ± 0.59 


0.38 ± 0.48 


-18.28 + 0.55 


62.0 ± 


5, 


.0 


69.0 




Blnpe 








17 


81377* 


149757 


C Oph 


16 


37 


09, 


.53 


-10 


34 


01.7 


7.12 ± 0.71 


13.07 ± 0.85 


25.44 + 0.72 


-9.0 ± 


5 


.5 


23.5 


348 


09.5Vnn 


15.9 


21.1 


0.16 


18 


82868 


152478 




16 


56 


08, 


.85 


-50 


40 


29.2 


4.34 ± 0.82 


-10.21 ± 0.84 


-9.55 + 0.62 


19.0 ± 


5 


.0 


30.3 




B3Vnpe 


7.9 


7.0 




19 


91599 


172488 




18 


40 


48, 


.06 


-08 


43 


07.5 


3.61 ± 1.16 


-9.64 ± 1.13 


-22.64 + 0.79 


34.1 ± 


1 


.2° 


44.7 




B0.5V 


12.7 


13.5 




20 


102274 


197911 




20 


43 


21, 


.62 


+63 


12 


32,9 


1.42 ± 0.62 


-13.72 ± 0.53 


-3.66 + 0.53 


-3.8 + 


5, 


.0 


46.1 




B5 








21 


109556* 


210839 


A Cep 


22 


11 


30, 


.58 


+59 


24 


52,3 


1.98 ± 0.46 


-7.22 ± 0.44 


-11.06 + 0.39 


-75.1 ± 


1, 


.2 


74.0 


214 


06I 


40.0 


64.6 


0.17' 


22 




J0826+2637 




8 


26 


51, 


.31 


+26 


37 


25,6 


2.6 


61 + 3 


-90 + 2 








P 


=0.53 


T = 4.92 


1.4" 






1 




,10835-4510 




8 


35 


20, 


.68 


-45 


10 


35.8 


2.0 


-48 ± 2 


35 + 1 










0.09 


0.01 


1.4" 






2 




J1115+5030 




11 


15 


38, 


.35 


+50 


30 


13.6 


1.9 


22 + 3 


-51 + 3 










1.65 


10.53 


1.4" 






4 




J1932+1059 




19 


32 


13, 


.87 


+ 10 


59 


31,8 


5.9 


99 + 6 


39 + 4 










0.22 


3.10 


1.4" 






8 




Geminga^ 




6 


33 


54, 


.15 


+ 17 


46 


12,9 


6.4 ± 1.7 


138 + 4 


97 + 4 














1.4" 






9 



Notes: a: Vrotsini from Slettebak, Wagner & Bertram (1997). b: Total mass for the binary: B5V (5.8 Mq) + F8V (1.1 M©). c: Vrad from Duflot et al. (1995). The 
Hipparcos Input Catalogue radial velocity is incorrect, d: Vrot sini from Bernacca & Perinotto (1970). e: We took the average Vrad from Bohannan & Garmany (1978), 
Garmany, Conti & Massey (1980), Stone (1982), and Gios & Bolton (1986). /: «rad, i-c, center of mass velocity, from Stickland et al. (1987). g: Wrotsini from Gies 
(1987). h: Spectral type of the secondary of t Ori from Stickland et al. i: Total mass for the binary t Ori. The individual masses are 22.9 AIq for the primary and 
14.9 Mq for the secondary, j: Wrot sin /' from Morse, Mathieu & Levine (1991). k: Wrotsini from Uesugi & Fukuda (1970). /: c from Kudritzki & Hummer (1990). m 
Assumed mass for main-sequence star instead of luminosity class IV. n: Viot sini from Brown & Verschueren (1997). o: Vrad from Gies & Bolton (1986). p: We take the 
characteristic mass for a neutron star, q: Position from Caraveo et al. (1998); parallax and proper motion from Caraveo et al. (1996). 
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3. A binary supernova in Upper Scorpius 

3.1. C, Oph and PSR J 1932+ 1059 

C Oph is a single 09.5Vnn star, and was first identified 
as a runaway originating in the Sco 0B2 association by 
Blaauw (19521;). Based on its proper motion, which points 
away from the association, its radial velocity, and the large 
space velocity (~30 km s^^), Blaauw suggested that C Oph 
might have formed in the center of the association ~3 Myr 
ago. Later investigations (e.g.. Paper I; Blaauw 1993; van 
Rensbergen et al. 1996D showed that ( Oph either became 
a runaway ~1 Myr ago in the Upper Scorpius subgroup of 
Sco 0B2, or 2-3 Myr ago in the U pper Centaurus Lupus 
subgroup (cf. de Zeeuw et al. 199£ ). 

If C Oph is a BSS runaway, as suggested by its high 
helium abundance (e = 0.16, corresponding to a mass 
fraction X = 0.577 of H) and large rotational velocity 
(348 km s^^), and if the binary dissociated after the su- 
pernova explosion, we might be able to identify the asso- 
ciated neutron star. None of the pulsars in Figure ^ was 
ever inside the Upper Centaurus Lupus subgroup, but two 
could have originated from the Upper Scorpius subgroup: 
PSR J1239+2453 and PSR J1932+1O590. 

We first consider PSR J1239-I-2453. Its estimated dis- 
tance is ~560 pc. It passed within about 20 pc of the 
Upper Scorpius region ~1 Myr ago if and only if its (un- 
known) radial velocity is large and positive (^650 km s~^). 
With a tangential velocity of ^300 km s~^ (the proper 
motion is 114 mas yr^^), the space velocity would have 
to be over 700 km s^^, which is uncomfortably large. 
Furthermore, while 1 Myr is consistent with the kine- 
matic age for ^ Oph, it is in conflict with the charac- 
teristic age (P/(2P) = 23 Myr) of the pulsar. The lat- 
ter is an uncertain age indicator, but the difference be- 
tween the two times is so large that we consider it un- 
likely that PSR J1239-f 2453 was associated with C Oph. 
The pulsar is currently at a Galactic latitude of 86°, i.e, at 
z ~560 pc above the Galactic plane. Typical z-oscillation 
periods of pulsars are of order 100 Myr (e.g., Blaauw & 
Ramachandran 1998), so that maximum height is reached 
after ~ 25 Myr. Taking the characteristic age at 

face value suggests the pulsar is near its maximum height 
above the plane, had a z-velocity of about 30 km s^^, and 
was not formed in the Upper Scorpius association (age 
^^5 Myr), but was born ~25 Myr ago in the Galactic plane 
outside the Solar neighbourhood. 

The path of the other pulsar, PSR J1932-f 1059 (ear- 
lier designation PSR B1929-fl0), also passed the Up- 
per Scorpius association some 1-2 Myr ago. The char- 
acteristic age of this pulsar is only ~3 Myr, consistent 
with the kinematic age of ^ Oph within the uncertain- 
ties. The present z-velocity of the pulsar (~40 km s^^ 



away from the Galactic plane) predicts a maximum dis- 
tance away from the plane of 680 pc and T1/4 ~ 28 Myr. 
The pulsar is presently located only ~10 pc below the 
plane. Since it presumably formed close to the plane, this 
means that PSR J1932-I-1059 either formed recently or 
well over 50 Myr ago. Considering that both the charac- 
teristic age and the typical pulsar ages (up to ^50 Myr) 
(Blaauw & Ramachandran 1998| ) are significantly smaller 
than ^50 Myr, we conclude that the pulsar formed re- 
cently. Upper Scorpius is the only site of star formation 
along the past trajectory of the pulsar. We thus consider 
PSR J1932-I-1059 a good candidate for the remnant of the 
supernova which caused the runaway nature of ^ Oph. 

3.2. Data 

Table ^ summarizes the data for ^ Oph and PSR 
J1932-f 1059. The radial velocity of the pulsar is unknown. 
The pulsar proper motion listed by Taylor et al. (1993 ) was 
calculated from timing measurements (Downs & Reichley 



Recently, Walter ( ^000| ) suggested that RX J185635-3754 
as another candidate neutron star that could have encountered 
^ Oph in the past. We show in Appendix B that this is unlikely. 



1983). More accurate proper motions can be obtained from 
VLBI observations; Campbell ( 1995| ) measured a provi- 
sional proper motion and parallax of PSR J1932-f 1059 
of {f^a*,fJ-s) — (96.7±1.7, 41.3±3.5) mas yr~^ and tt = 
5±1.5 mas, respectively, including a full covariance ma- 
trix. These measurements are in good agreement with 
those of Taylor et al. ( |1993| ; see Table | and Figure |). 

3.3. Simulations 

Our hypothesis is that C Oph and PSR J1932-I-1059 are 
the remains of a binary system in Upper Scorpius which 
became unbound when one of the components exploded 
as a supernova. Support for this hypothesis would be to 
find both objects at the same position at the same time 
in the past. Our approach is to calculate their past orbits 
and simultaneously determine the separation between the 
two objects, -Dmin(''"), as a function of time, r. We define 
-Dniin(''') as I a; (J oph — a^puisarl, where Xj is the position 
of object j. We consider the time tq at which Z3min(''') 
reaches a minimum to be the kinematic age. To take the 
errors in the observables into account we calculate a large 
set of orbits, sampling the parameter space defined by 
the errors. We use the Taylor et al. ( |1993| ) proper mo- 
tion for the pulsar. The errors in the positions of the 
runaway and the pulsar are negligible, and those in the 
proper motions of the two objects and in the parallax of 
the runaway are modest (< 10%). However, the radial- 
velocity error of ( Oph is considerable (5 km s~^). The 
distance to the pulsar has a significant error, and its ra- 
dial velocity is unknown. Accordingly, we first determine 
the region in the (TTpuisan ■^rad, pulsar) parameter space for 
which the pulsar approaches the runaway when we re- 
trace both orbits. Sampling a grid in (TTpuisar, Wrad, pulsar) 
while keeping the other parameters fixed, we find that, for 

2 ^ TTpulsar ^ 6 maS and 100 ^ Wrad^pulsar ^ 300 km s~\ 
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Fig. 3. Left: Distribution of minimum separations, 
Anin(To), between C Oph and PSR J1932-I-1059. The solid 
line denotes the expected distribution of i^min, see §3.3. 
Right: Distribution of the times tq at which the minimum 
separation was reached. 



the motions of the pulsar and the runaway are such that 
their separation decreases as one goes back in time. 

Adopting TTpuisar = 4±2 mas and Wrad, pulsar = 200±50 
km s~^, we calculate three million orbits for the pulsar and 
the runaway. Considering that the pulsar proper-motion 
errors might be underestimated (Campbell et al. 1996; 
Hartman 1997), we increased them by a factor of two. For 
each run we create a set of positions and velocities for 
the runaway and the pulsar consistent with the (modified 
3a) errors on the observables. We also calculated the or- 
bit of the Upper Scorpius association back in time, using 
the mean position and velocity derived by de Zeeuw et al. 
( |1999| , their table 2). 30 822 of these simulations resulted 
in a minimum separation between the pulsar and the run- 
away of less than 10 pc. In 4 214 simulations the pulsar and 
the runaway had a minimum separation less than 10 pc 
and were both situated within 10 pc of the center of the 
association (the smallest minimum separation found was 
0.35 pc). Thus, only a small fraction (0.14%) of the sim- 
ulations is consistent with the hypothesis that the pulsar 
and the runaway were once, ~1 Myr ago, close together 
within the Upper Scorpius association. We now show that 
given the measurement uncertainties, this low fraction is 
perfectly consistent with the two objects being in one lo- 
cation in the past. 

Figure |^ shows the distribution of the minimum sepa- 
rations, Djnin{To), and the kinematic ages, tq, of the 4 214 
simulations mentioned above. The lack of simulations 
which yield very small minimum absolute separations is 
due to the three-dimensional nature of the problem. Con- 
sider the following case: two objects are located at exactly 
the same position in space, e.g., the binary containing the 
pulsar progenitor and the runaway. However, the position 
measurement of each of these objects has an associated 
typical error. The distribution of the absolute separation 
between the objects, obtained from repeated measure- 
ments of the positions of both objects, can be calculated 
analytically for a Gaussian distribution of errors (see 
Appendix A). The solid line in Figure || shows the result 



for an adopted distance measurement error of 2.5 pc, 
which agrees very well with our simulations. The peculiar 
statistical properties of the sample of successful simula- 
tions make it difhcult to give a simple argument to derive 
the value of 2.5 pc from the uncertainties in the kinematic 
properties of the runaway and the pulsar. We suspect that 
the disagreement between the solid line and histogram for 
separations >6 pc is most likely due to a slight mismatch 
between the model and the actual situation. Even so. 
Figure ^ shows that due to measurement errors, very 
few simulations will produce a small observed minimum 
separation, even when the intrinsic separation is zero. 

Figure ^ shows the astrometric parameters of the pul- 
sar and the runaway at the start of the orbit integration, 
i.e., the 'present' observables, for the simulations which 
result in a minimum separation less than 10 pc occurring 
within Upper Scorpius. The parameters of C Oph show no 
correlations except between the parallax and the radial 
velocity. This is expected due to the degeneracy of these 
two quantities (a change in stellar distance, depending on 
whether it increases or decreases the separation between 
the star and the association, can be compensated for by a 
larger or smaller radial velocity, respectively). 

The parameters of the pulsar behave very differently. 
In addition to the tt vs. Wrad correlation, we find that the 
parallax is also correlated with both of the proper mo- 
tion components. As a result, the proper motion compo- 
nents are correlated with each other. This means that only 
a subset of the full parameter space defined by the six- 
dimensional error ellipsoid of the pulsar fulfills the require- 
ment that the pulsar and the runaway meet. Furthermore, 
if they met, then we know the radial velocity of the pulsar 
for each assumed value of its distance. A reliable distance 
determination would thus yield an astrometric radial ve- 
locity. The current best distance estimate of the pulsar 
deri ved fr om VLBI measurements, tt — 5±1.5 mas (Camp- 
bell 1995), predicts a radial velocity of 100-200 km s~^. 
This radial velocity is comparable to the tangential veloc- 
ity: ^100 km s^^ (for tt — 5 mas). 

The pulsar proper motions of the 4 214 successful sim- 
ulations are shown in Figure |^, together with the proper- 
motion measurements of Lyne, Anderson & Salter (1982) 
(dot-dash line), Taylor et al. ( |l993| ) (solid line), and 
Campbell ( 1995 ) (dashed line). The measurements show 
a reasonable spread, reflecting the difficulty in obtaining 
pulsar proper motions, but are consistent, within 3cr, with 
the proper motions predicted by the simulations. 

3.4- Interpretation 

The observed astrometric and spectroscopic parameters 
of C Oph and PSR J1932-I-1059 are consistent with the 
assumption that these objects were very close together 
~1 Myr ago (Figure |^). At that time both were within 
the boundary of Upper Scorpius (Figure ^ , which has a 
nuclear age of ^--^5 Myr (de Geus, de Zeeuw & Lub 1989). 
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Fig. 4. Astrometric parameters of the runaway ( Oph and the pulsar PSR J1932+1059 at the start of each of the 
4 214 simulations for which the minimum separation between the runaway and the pulsar was less than 10 pc, and 
both of them were within 10 pc of the center of Upper Scorpius sometime in the past. 



Several characteristics of ( Oph and Upper Scorpius 
support the interpretation that the runaway and the pul- 
sar must have been produced by the BSS. (i) The HI dis- 
tribution in the direction of Upper Scorpius shows an ex- 



panding shell-like structure. De Geus (1992) argued that 
the energy output of the stellar winds of the massive stars 
in Upper Scorpius is two orders of magnitudes too low to 
account for the kinematics of the HI shell, and proposed 
that a supernova explosion created the Upper Scorpius HI 
shell, (ii) Based on the present-day mass function, de Geus 
showed that the initial population of Upper Scorpius most 
likely contained one star or binary more than the present 
population. The estimated mass of this additional object is 
~40 Mq. Stars of this mass have main-sequence lifetimes 
of ~4 Myr, and end their lives in a supernova explosion. 
Since the association has an age of ^^5 Myr, the supernova 
explosion might have taken place ~1 Myr ago. (iii) The 
characteristics of C Oph are indicative of close binary evo- 
lution. The helium abundance is large, and the star has 
a large rotational velocity (see These facts make it 
very likely that the same supernova event in Upper Scor- 
pius created PSR J1932-I-1059 and endowed ( Oph with 
its large velocity. 

It could be that C Oph and PSR J1932-t-1059 are not 
related. This would imply that the pulsar originated in 
Upper Scorpius ^1 Myr ago and that C, Oph obtained 
its large velocity either in a separate BSS event in Upper 
Scorpius, or in Upper Centaurus Lupus, '^S Myr ago, as 




60 80 100 120 

f^a*, pulsar |nias yr-' j 

Fig. 5. Proper motions of the pulsar PSR J1932-K1059 
at th e start of the 4 214 successful simulations (dots; see 
§3.3), and the proper-motion measurements of the pulsar: 
dot-dash line denotes Lyne, Anderson & Salter (1982), 
solid line denotes Taylor, Manchester & Lyne (1993), and 
the dashed line denotes Campbell (1995). The contours 
indicate the 1, 2, and 3(t confidence levels. 



suggested by van Rensbergen et al. (1996). In the latter 



case, it would also have to be formed by the BSS, because 
of its high helium abundance, large rotational velocity, 
and the ~10 Myr difference between the age of Upper 
Centaurus Lupus and the kinematic age of C Oph (§|l]). 
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The simulations described in §B.3 provide the veloci- 
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Fig. 6. The orbits of C Oph, PSR J1932+1059, and Up- 
per Scorpius. The present positions are denoted by a star 
for the runaway, a filled circle for the pulsar, and by a 
filled square for the association. The top panel shows the 
distance vs. Galactic longitude of the stars. The bottom 
panel shows the orbits projected on the sky in Galactic 
coordinates. The small open circles in the bottom panel 
denote the present-day positions of the O, B, and A- type 
members of Upper Scorpius, taken from de Zeeuw et al. 
( |1999D . The large circle denotes the position of the asso- 
ciation at the time of the supernova explosion, and has a 
10 pc radius. This figure assumes a set of space motions 
consistent with the common origin hypothesis. 



Given the small probability of finding a runaway star and 
a pulsar with orbits that cross, and with both objects at 
the point of intersection at the same time, we conclude 
that C Oph and PSR J1932-H1059 were once part of the 
same close binary in Upper Scorpius, providing the first 
direct evidence for the generation of a single runaway star 
by the BSS. 

3.5. Pulsar kick velocity 

If C Oph and PSR J1932-hl059 were once part of a bi- 
nary, then we can derive a number of properties of this 
system. For example, the true age of the pulsar must be 
the kinematic age of 1 Myr (as compared to the character- 
istic age estimate of 3 Myr). It follows that, if no glitches 
occurred, the pulsar had a period of 0.18 seconds at birth, 
as compared to the current period of 0.22 seconds. 

The velocity distribution of the pulsar population is 
much broader (a fewx 100 km s^^) than that of the pulsar 
progenitors (a fewx 10 km s^^). The mechanism responsi- 
ble for this additional velocity (the 'kick velocity' Wkick) is 



not well understood (e.g., Lai 1999). The kick velocity is 



ties of the runaway star, the pulsar, and the association at 
the present time. This makes it possible to determine i»kick 
of the neutron star. For the 4 214 successful runs we find 
that the average velocities with respect to Upper Scorpius 
are: t;c Oph = (-6.4±4.2, 33.8±1.4, 5.8±2.0) km s^^ 
and Vpuisar = (48.6±21.7,222.9±36.1,-70.7±8.4) 
km s""'^ in Galactic Cartesian coordinates {U, V, W)^ 

To derive t^kick we consider a binary with components 
of mass Ml and M2 in a circular orbit, in which the first 
component (stari) explodes as a supernova and creates 
a neutron star. At the time of the explosion, stari is 
the least massive component of the binary, due to the 
prior mass transfer phase, and is most likely a helium 
star. The rapidly expanding supernova shell, with mass 
AM = Ml- Mn, where Mn = 1.4 Mq is the typical mass 
a neutron star, will quickly leave the binary system. The 
shell has a net velocity equal to the orbital velocity of 
stari at the moment of the explosion (vi). A net amount of 
momentum (AM x vi) is thus extracted from the system 
and the binary reacts by moving in the opposite direction 
with a velocity v = ~{AM x wi)/(M2-|-M„), the so-called 
'recoil velocity'. The binary will remain bound after the 
explosion because less than half of the total mass of the 
system is expelled (Mi < M2; cf. Paper I). However, if the 
neutron star receives a kick in the supernova explosion the 
binary might dissociate, depending on the direction and 
magnitude of the kick velocity. We simulate this by using 
a simple orbit integrator for two bodies. We determine 
the semi-major axis and orbital velocities assuming the 
binary has a circular orbit and masses Mi = 5 Mq and 
M2 = 15 Mq {( Oph). We then change the mass of stari to 
M„ and add a kick- velocity to its orbital velocity. We start 
the integration at this point and try to reproduce the ob- 
served velocity of ( Oph, the pulsar, and the angle between 
the two velocity vectors (35°). It turns out that a kick ve- 
locity of order 350 km s~^ is needed in a direction almost 
opposite to the orbital velocity of stari 's prior to the ex- 
plosion. This value is in good agreement with the average 
pulsar kick velocity found by Hartman (1997) and Hansen 
& Phinney ( |l997| ). The current velocity of the pulsar, 
~240 km s^^, is more than 100 km s^^ smaller than the 
kick it acquired. Our simulations show that this decelera- 
tion is due to the gravitational pull of Oph on the pulsar. 

The mass of ( Oph used in the above estimate is con- 



sistent with the calibration of Schmidt-Kaler (1982). The 



more recent mass calibration of Vanbeveren, Van Rensber 
gen & De Loore ( |l998| ) suggests 21 Mq (Table |) 



would increase the inferred kick velocity to 



4. A dynamical ejection in Orion 



This 
400 km s-i. 



most likely due to asymmetries in the core of a star just 
before, or during, the supernova explosion. 



^ U points in the direction of the Galactic centre, V points 
in the direction of Galactic rotation, and W points towards the 
North Galactic pole. 
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Fig. 7. Contours of minimum separation for the pairs 
AE AuT-fi Col {top), AE Aur-6 Ori {middle), and ^ Col- 
L Ori {bottom). The contours are spaced every 1 pc with 
the outermost contour being 10 pc. The ordinates repre- 
sent straight lines in the distance-distance plane, defined 
as indicated in the top right of each panel. 

AE Aurigae & /i Columbae 



Blaauw & Morgan (1954) drew attention to the isolated 
stars AE Aur (09.5V) and ^^ Col (O9.5V/B0V), which 
move away from the Orion star-forming region (e.g., Mc- 



Caughrean & Burkert 2000) in almost opposite directions 
with comparable space velocities of ~100 km s^^ (Fig- 
ure |l^, stars 5 and 6 in Table ^ . Blaauw & Morgan sug- 
gested that . . the stars were formed in the same physical 
process 2.6 million years ago and that this took place in 
the neighborhood of the Orion Nebula. ' The past orbits of 
AE Aur and /i Col intersect on the sky near the location 
of the massive highly-eccentric double-lined s pectro scopic 
binary l Ori (09III +B1I II, see Stickland et al. |1987D . This 
led Gies & Bolton ( 1986 ) to suggest that the two runaways 
resulted from a dynamical interaction also involving t Ori: 
i Ori is the surviving binary of a binary-binary colli- 
sion that ejected both AE Aur and fi Col. ' 

4.2. Data 

Table ^ hsts the data for AE Aur, Col, and l Ori. We 
adopt Stickland's et al. ( |1987| ) radial velocity for t Ori 
{v^). For the radial- velocity errors for AE Aur and /i Col 



we use the largest errors quoted in either the Catalogue 
de Vitesses Radiales Moyennes Stellaires (Barbier-Brossat 
1989), the Hipparcos Input Catalogue (Turon et al. 1992), 
the Wilson-Evans-Batten Catalogue (Duflot, Figon & 
Meyssonnier |1995|), or in the SIMBAD database. 



4.3. Simulations 

To investigate the hypothesis that the three stellar sys- 
tems, AE Aur, fi Col, and t Ori, were involved in a binary- 
binary encounter, we retrace their orbits back in time to 
find the minimum separation between them. As in we 
explore the parameter space determined by the errors of, 
and correlations between, the observables. 

Even with the unprecedented accuracy in trigono- 
metric parallaxes obtained by the Hipparcos satellite, 
the errors on the individual distances are rather large: 
Dae Aur = UGtlTi PC, coi = 597^^11 pc, D, ori = 
406l^gg pc. We therefore first determine which distances 
are most likely to agree with our hypothesis, and then 
study the effect of the measurement errors on the other 
observables. For each pair of stars. Figure shows contours 
of minimum separation between the respective orbits as a 
function of their present distances. The distances of the 
stars for which the orbits have a small minimum separa- 
tion are strongly correlated, i.e., if the distance of star i 
increases that of star j also needs to increase to obtain a 
small minimum separation. We therefore choose to show 
the contours of constant minimum separation with respect 
to this correlation. The vertical axes thus show offsets from 
the straight line in the distance vs. distance plane defined 
by the equation in the top right of each panel. 

We start each simulation with a set of positions and ve- 
locities which are in agreement with the observed param- 
eters and their covariance matrices (<3f7). Furthermore, 
we require the distances of the stars to fall within the 
10 pc minimum-separation contours of Figure 0. We then 
calculate the orbits of AE Aur, n Col, and l Ori. We define 
the separation between the three stellar systems, £'min(T), 
as the maximum deviation of the objects from their aver- 
age position, i.e., Dnun{T) = maxjccj — x\ for j = AE Aur, 
fj. Col, and l Ori, where x = ^{xae Am +a^p Coi + a^t Ori) is 
the mean position and Xj the position of star j. The time 
To at which I?rnin(''') rcaches a minimum is considered to 
be the time of the encounter, i.e., the kinematic age. 

We computed 2.5 million orbits, of which 114 yielded 
Dnun{To) < 1 pc with tq = 2-3 Myr. One of the simula- 
tions resulted in a minimum separation of 0.019 pc which 
is equal to 4 000 AU (Figure ^) . The small number of sim- 
ulations with small minimum separations is due to (i) the 
large number of parameters involved (i.e., 18) and (ii) the 
three-dimensional nature of the problem (cf. §3.3). 

We have numerically determined the distribution of 
the minimum separations Dmin of three points drawn from 
a three-dimensional Gaussian error distribution (the ana- 
lytic results of Appendix A are valid only for two Gaus- 
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identical to the mean velocity Vdus of the parent cluster. 
Then 
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Fig. 8. Left: Distribution of minimum separations be- 
tween AE Aur, fi Col, and l Ori, Dmin(T'o), of 10 000 Monte 
Carlo simulations of the stellar encounter. Right: The 
-Dniin(To) distribution for three randomly drawn points 
from three spherical Gaussians with standard deviations 
of cr = 4 pc (solid line and shaded) , cr = 2 pc (dotted hue) , 
and CT = 6 pc (dashed line). The -Dmin(''"o) distribution for 
cr = 4 pc is a good representation of the distribution in 
the left panel. Four pc is the typical spread in the end po- 
sitions of the orbits due to the errors on the present day 
velocity (~2 km s~^). The dotted and dashed histograms 
have been normalized such that their shapes can be com- 
pared with the solid histogram. 



sians). We randomly draw three points from three spheri- 
cal three-dimensional Gaussians (with standard deviation 
a) and determine Umin- The Gaussians have the same 
mean positions. The resulting distribution for cr = 4 pc 
resembles the real one remarkably well (Figure ||). A 
distance uncertainty of four pc is consistent with the 
~2 km s~^ uncertainties in the velocities of the runaways 
and t Ori: 2 km over ~2 Myr results in a displacement 
of ~4 pc. Thus, the data and their errors are consistent 
with the hypothesis that ~2.5 Myr ago AE Aur, fi Col, 
and i Ori were in the same small region of space. 

4-. 4- Interpretation 

The nominal observed properties of the runaway stars 
AE Aur and /i Col and the binary l Ori are consistent with 
a common origin ^2.5 Myr ago. The most likely mecha- 
nism that created the large velocities of the runaways and 
the high eccentricity of the t Ori binary is a binary-binary 
encounter, as suggested by Gies & Bolton ( |1986| ). The nor- 
mal rotational velocities of both runaways (25 km s~^ and 
111 km s~^) and the normal helium abundance of AE Aur 
(Table H see Blaauw [l993[ figure 6) also suggest that these 
runaways were formed by the dynamical ejection scenario. 
The helium abundance of fi Col is unknown. 

4-. 5. Parent cluster 

To find the cluster, or region of space, where the encounter 
between AE Aur, n Col, and t Ori took place we assume 
that the center of mass velocity of the three objects is 



'^clus 



(2) 



for j = AE Aur, /i Col, and t Ori. For each star we es- 
timate the mass by interpolating the mass vs. spectral- 
type calibration of Schmidt-Kaler (1982, table 23). We 
obtain 15.9 Mq for AE Aur and ^ Col, and 22.9 Mq and 
14.9 Mq for the primary and secondary of i Ori, respec- 
tively (37.8 Mq for the binary system). We use the cluster 
velocity (f cius) and the mean position of the three stellar 
systems at the moment of the encounter to integrate the 
orbit of the ensemble of stars tq Myr into the future. The 
position and velocity at the end of this integration should 
coincide with the present-day properties of the parent clus- 
ter. We extend the Monte Carlo simulations described in 



§4.3 to include the integration of the orbit of the 'cluster' 
forward in time. Figures I and[l^ summarize the results. 
Panel b of Figure |l^ shows that the distances of the three 
stars and the predicted cluster distance are tightly corre- 
lated (see also Figure 0); all distances increase when the 
cluster distance increases. A consequence of this tight cor- 
relation is that as soon as the distance to one of the objects 
is known, all other distances are fixed. 



4.5.1. Biases and measurement errors 

Two effects influence the mean cluster properties as pre- 
dicted by the Monte Carlo simulations. First, it is eas- 
ier to hit a target from close by than from far away, i.e., 
a larger range of velocities (within the errors) is consis- 
tent with the encounter hypothesis when the distance be- 
tween the star and the encounter point is small (the 'aim- 
ing effect'). We simulate this effect in the following way. 
We assume a range of cluster distances, 350-500 pc. For 
each distance we use Figure ^ (the gray dots in the first 
row) to determine the other phase-space coordinates of 
the parent (position on the sky, proper motion, and ra- 
dial velocity). With these 'observables' we calculate the 
three-dimensional velocity of the cluster, corrected for So- 
lar motion, and determine its position at a time tq (see 
first row in Figure ||) in the past. We neglect the vari- 
ation of the Galactic potential, ignore Galactic rotation, 
and use the linear velocity, to speed up the calculations. 
This past position of the cluster combined with the present 
three-dimensional positions of AE Aur, fi Col, and l Ori 
(based on the present positions on the sky and the dis- 
tances from Figure |l^ panel b) gives the velocities of the 
three stellar systems today, using tq as the time difference. 
These 'observed' properties are then used as input for the 
Monte Carlo simulations described above to investigate 
the influence of the aiming effect on the predicted cluster 
distance. The circles in Figure |l^ panel c display the bias 
in the cluster distance. 
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Fig. 9. Properties of the parent cluster of the runaways AE Aur and /i Col and the binary l Ori obtained from our 
Monte Carlo simulations. First row: Cluster properties plotted vs. the cluster distance. The grey dots denote the 
median values of the cluster properties for distance-bins of 25 pc. Second row: Histograms of the predicted cluster 
properties. The tick marks on the vertical axis have a spacing of 1000. Third row: Histograms of the predicted cluster 
properties when the mass of Col is changed by —1 Mq. Note that the distance, time, and radial- velocity histograms 
do not change significantly. The tick-mark spacing along the vertical axis is similar to that in the second row. 



Secondly, the trigonometric distance of /i Col, 



Col — 397^^7° pc, is smaller {2a) th an th e observed 



photometric distance, ^--^750 pc (e.g., Gies 1987). The pho- 
tometric distance is reliable since ^ Col is located in a 
region free of interstellar absorption. This difference be- 
tween the trigonometric distance and the 'real' distance 
results in an additional bias towards smaller distances for 
the stars and the cluster. In our Monte Carlo simulation 
we draw the parallaxes, like all other observables, from a 
Gaussian centred on the observed value and with a width 
equal to the observed error. For the Hipparcos distance 
of /i Col this means that less than ^10% of the random 
realizations will be consistent with the photometric dis- 
tance]^. And because the distances of the three stellar sys- 
tems and the cluster are correlated (see Figure |l^ panel 
6), the other stars also need to be at smaller distances 
for the encounter to take place. This effect will result in a 
mean cluster distance (the mean of the Monte Carlo simu- 
lations) which is underestimated. We simulated this effect 
in a similar manner as the aiming effect. The results on 
the mean cluster distance in the Monte Carlo simulations, 



® The Hipparcos parallax of /i Col (2.52±0.55 mas) deviates 
more than 1.5(t from the photometric parallax (~1.3 mas). 
The random realisations of vr^ coi which are consistent with 
TTphot thus need to fall outside the — 1.5cr confidence level which 
occurs in fewer than 10% of the cases. 



aiming effect and the parallax of /i Col, are shown as the 
triangles in Figure M panel c. 



4.5.2. Cluster properties and identification 

Taking the biases on the cluster distance into account, we 
reconstruct the present-day properties of the parent clus- 
ter of the stars AE Aur, Col, and l Ori. The mean clus- 
ter distance from our Monte Carlo simulations is 339 pc 
(right most panel in the second row of Figure |^). The clus- 
ter distance corrected for biases is 425-450 pc. Using this 
distance we determine the other properties of the cluster 
(first row of Figure ^), and summarize them in Table ^ 
Figure ^ (right most panel of the first row) indicates that 
the encounter happened 2.5 Myr ago; this obviously is a 
lower limit to the age of the cluster. Figure |l^ shows the 
region of the sky where the parent cluster should be lo- 
cated: the Orion Nebula. The black contours in the bottom 
panel show the distribution of (^, h) of the parent cluster 
obtained from the Monte Carlo simulations. 

Of all clusters in this active star-forming region the 
Trapezium cluster (NGC 1976) is the most likely parent 
cluster for the following reasons. 



1: The cluster is young. Palla & Stabler ( |1999| ) find a 
mean age of ~2 Myr based on theoretical pre-main- 
sequence tracks, and established that the first stars 
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Fig. 10. Properties of the parent cluster of the runaways AE Aur and /i Col and the binary t Ori. Panel a: Proper 
motions and their errors (grey circles) for all stars in the Tycho 2 Catalogue (filled triangles) within an area of 0?4 by 
0?4 centred on the Trapezium cluster. The large grey dot denotes the average of the predicted cluster proper motion 
for the Monte Carlo simulations. Panel b: Distances of the runaway stars as a function of the cluster distance in the 
Monte Carlo simulations. The different grey scales are labeled in the panel. The filled circles and their error bars 
denote the observed distances of the stars derived from the Hipparcos parallaxes (prlx) and the open circles denote 
the distances derived from photometry (phot) (Gies 1987). Panel c: the biases on the predicted cluster distance as 
discussed in §4.5.1, The filled and open symbols denote the mean and median, respectively, of the cluster distance 
distributions. The circles include only the first bias, the aiming effect, and the triangles include both the aiming 
effect and the 'incorrect' Hipparcos parallax of /i Col. For clarity, the circles and triangles are displaced —3 and 3 pc, 
respectively. The dotted line indicates the mean cluster distance based on the Monte Carlo simulation. 



Table 4. Predicted properties of the parent cluster of the 
stellar systems AE Aur, /x Col, and i Ori: the cluster dis- 
tance -Dparont, the sky position in equatorial {a,S) and 
Galactic {6,b) coordinates, the proper motions (/iQ*,/!^) 
and {fJ.it: , Hb) , and the radial velocity frad- The predicted 
distances of the runaways and t Ori if there was an 
encounter: Dae Aur = 430 pc, coi — 600 pc, and 
D, Ori = 440 pc. 

Col - 1 Mq 



-Dparciit 

{a, 5} 



425-450 
(84?0,-5?8) 
(1.7,-0.8) 
(209?4,-19?4) 
(1.3,1.2) 
28.3 



425-450 
(83°9,-5?2) 
(1.7,-0.2) 
(208?9,-19?2) 
(0.9,1.4) 
27.6 



mas yr 



mas yr 
km s^^ 



formed not more than 5 Myr ago. Thus, the Trapezium 
is old enough to have produced the runaways. 
The Trapezium is one of the most massive, dense clus- 
ters in the Solar neighbourhood. Estimates for the stel- 
lar density are > 20 000 stars pc^^ for the inner 0.1-0.3 
pc (e.g., McCaughrean & Stauffer 1994; Hillenbrand & 



Hartmann 1998). These high stellar densities favor dy- 
namical interactions within the cluster core. 
The Trapezium shows a strong mass segregation (Zin- 
necker, McCaughrean & Wilking 1993; Hillenbrand & 



Hartmann |1998|) 
than 



Five of the six stars more massive 
are in the centre. This concentration of 



IOMq 

massive stars increases the probability for dynamical 
interactions between these stars. 



The binary fraction in the Trapezium cluster is at least 
as high as that of the Solar-type field stars, i.e., ^60% 
(Prosser etal. |l994t Petr e t al. [1998 ^ Simon, Cose 



Beck |1999| ; Weigelt et al. |1999| ). This means that 



enough binary systems are available for binary-binary 
or binary-single-star interactions to become efficient in 
expelling stars from the cluster. 

The mean astrometric properties and the radial ve- 
locity of the Trapezium agree perfectly with those pre- 
dicted by our Monte Carlo simulation. The distance to the 



Trapezium is estimated to be 450-500 pc (Walker 1969; 
Warren & Hesser |T977a| p^77b| , |l97|; Genzel & Stutzki 



1989); we predict 425-450 pc. The observed radial velocity 
of the Trapezium is 23-25 km s^^ (Johnson p_965 ; Warren 



& Hesser pW7a|, pW7b| ; Abt, Wang & Cardona ]l99lt Mor- 
rell & Levato 1991); we predict ~28 km s""'^. The absolute 
proper motion of the Trapezium is ill-det ermin ed, but is 
known to be small (e.g., de Zeeuw et al. 1999). We col- 



lected all stars, within a 0?4 by 0?4 region centred on th e 
Trapezium, based on the Tycho 2 Catalogue (H0g ^000 ), 
and plot the proper motions in Figure |lO| a. The proper 
motions agree with the predicted cluster proper motion. 

Table ^ shows that the predicted position on the sky 
of the parent cluster does not fully agree with the posi- 
tion of the Trapezium (see Figure |l^) . Here it is important 
to remember that we did not allow for any errors on the 
stellar masses used in eq. (||). We investigate the effect of 
mass errors by changing the masses and running a new set 
of Monte Carlo simulations. We find that (i) the results 
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are insensitive to the mass of l Ori: a change as large as 
±5 Mq produces no noticeable change in the cluster prop- 
erties, and (ii) the sky position of the parent cluster and 
its proper motion depend on the mass ratio of AE Aur 
and /z Col. Changing the mass of Col by —1 Mq or 
the mass of AE Aur by +1 Mq shifts the predicted sky 
position of the parent cluster to that of the Trapezium 
cluster (Figure |Tl|) . A mass change in the other direction, 
+1 Mq for Col and —1 Mq for AE Aur, creates a similar 
shift in the opposite direction. There are indications from 
spectral-type determinations that ji Col is indeed slightly 
less massive than AE Aur. Most spectral- type determina- 
tions of Col give 09.5V; however, Blaau w fc Morgan 
( |1954D and Paper I quote BOV and Houk ( [l982[ ) quotes 
BlIV/V. 
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Fig. 11. Top & middle: Orbits, calculated back in time, of the runaways AE Aur (dotted line) and /i Col (solid line) 
and the binary l Ori for one of the Monte Carlo simulations described in the text. The top panel shows the distance vs. 
Galactic longitude of the stars. The middle panel shows the orbits projected on the sky in Galactic coordinates. The 
starred symbols depict the present position of the three stars. The stars met ^2.5 Myr ago. Using conservation of linear 
momentum, the orbit of the parent cluster (grey solid line, see blow up) is calculated from the time of the assumed 
encounter to the present. The large circles denote all stars in the Hipparcos Catalogue brighter than V = 3.5 mag; 
filled circles denote O and B stars, open circles denote stars of other spectral type. The small circles denote the O 
and B type stars with 3.5 mag < V < 5 mag (cf. figure 1 in Blaauw & Morgan 1954). The Orion constellation is 
indicated for reference. Bottom.: The predicted position of the parent cluster (black contours) together with all stars in 
the Tycho Catalogue (ESA 1997) in the field down to V = 12.4 mag. The size of the symbols scales with magnitude; 
the brightest star is l Ori. The Trapezium and l Ori are indicated. The black and dark grey lines are the past orbits of 
L Ori and the Trapezium, respectively (see top panel). The triangle denotes the predicted present-day position of the 
parent cluster for this particular simulation. The grey contours display the IRAS 100 micron flux map, and mainly 
outline the Orion Nebula. 
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We note that the caUbration of Vanbeveren et al. 
( |1998D gives a mass of 38.6 M© for t Ori, similar to that 
found with the Schmidt-Kaler calibration, but increases 
the masses of AE Aur and /i Col to 21.1 Mq. This does 
not change our results, as it is the ratio of the runaway 
masses that determines the predicted current position of 
the parent cluster. 

In summary, the position, distance, proper motion, and 
radial velocity of the Trapezium cluster fall within the 
range predicted by our Monte Carlo simulations. Further- 
more, the youth, extreme stellar density, mass segregation, 
and the high binary fraction make it the best candidate 
for the parent cluster of the runaways AE Aur and /i Col 
and the binary t Ori. Finally, it is the only likely candidate 
in this region of the sky. 

5. 53 Ari, ^ Per, C. Pup, and A Cep 

The previous sections gave a specific example of each for- 
mation mechanism for runaway stars, and described our 
orbit retracing methods in detail. We now consider the 
three other classical runaways, as well as C Pup, and we 
discuss the likely formation mechanisms. The results are 
summarized in Table ^ 



450 



5.1. 53 Arietis & Orion OBI 



(star 2) 



Blaauw (1956) classified 53 Arietis (HIP 14514) as a run- 
away star based on its proper motion which is directed 
away from the Orion association. He deduced that 53 Ari 
left the Orion association ~4.8 Myr ago and predicted the 
radial velocity of the star (at that time unknown) to be 
~18 km s^^. The discrepancy in kinematic ages of AE Aur 
and /i Col (§||) and that of 53 Ari indicates that it is not 
related to the same event that created AE Aur and /i Col 
(sjj) but is another runaway from the Orion star-forming 
region. Table |^ lists the observables of 53 Ari[|. 

The Ori OBI association has four subgroups: a, b, c, 



and d (Blaauw 1964 ; Brown, de Geus & de Zeeuw 1994). 
We do not consider subgroup d (the Trapezium) as a pos- 
sible parent group of 53 Ari, since this subgroup is younger 
than the runaway {^^. The ages of the other subgroups 
are: 8-12 Myr for subgroup a, 2-5 Myr for subgroup b, and 
~4 Myr for subgroup c (Warren & Hesser 1977a, 1977b; 
Brown et al 



1994) 



Simulations: We performed a set of simulations as in §g, 
retracing orbits for each subgroup (a, b, c) . The kinematic 

The Hipparcos Input Catalogue (Turon et al. 1992) lists 
the radial velocity from the Catalog ue de Vitesses Radiales 
Moyennes Stellaires (Barbier-Brossat 1989) which is incorrect 
(—8.5 km s~^). The radial velocity in the Troisieme Cata- 
logue Bibliographique de Vi tesses Radiales Stellaires (Barbier- 
Brossat, Petit & Figon 1994 ) is also incorrect. They list a radial 
velocity correc ted for Solar motion of 15.3 km s~^ adopted from 
Sterken (1988). The uncorrected radial velocity is 24.2 km s~^. 
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Fig. 12. Predicted distance (left) and kinematic age 
(right) of 53 Arietis as a function of the distance of the 
parent association: Ori OBI subgroup a (top), Ori OBI 
subgroup b (middle), and Ori OBI subgroup c (bottom). 
The Hipparcos distance and its la error are indicated in 
the left panels. Distance estimated by Warren & Hesser 
(1977a; 1977b, WH77), Brown et al. (1994, B94), and de 
Zeeuw et al. (1999, Z99) are indicated in the middle pan- 
els. 



age of 53 Ari from subgroup a is ~4.3 Myr (Figure |l2|). 
This means that the subgroup was ~6 Myr old when 
53 Ari became a runaway star. This very likely rules out 
the DES as the formation mechanism (see [Q). However, 
there is little direct evidence in favor of the BSS. The 
helium abundance of 53 Ari is unknown and its observed 
rotational velocity is small (v^ot sin i = 10 km s-i), but 
this could be caused by a near pole-on orientation. We 
did not find a neutron star associated with 53 Ari, but 
our samp ling of the nearby compact objects is severely 
limited (§|t|). 

If subgroup b is the parent association the kinematic 
age for 53 Ari is ^-^4.8 Myr. This is comparable to the 
canonical age of the subgroup, and excludes the BSS as 
a production mechanism for 53 Ari (see §|l|). If Ori OBI 
b is the parent group of 53 Ari then the kinematic age is 
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~4.8 Myr and the formation mechanism is most hkely the 
DES. However, the most recent age determination (Brown 
et al. |1994| ) gives 1.7 ± 1.1 Myr. If Ori OBI b is indeed 
this young then the subgroup is younger than 53 Ari and 
cannot be the parent group. 

For subgroup c we find that the minimum separation 
between the subgroup centre and the runaway was never 
smaller than 15 pc, while the simulations for the other two 
subgroups a and b yield minimum separations as small as 
1 pc. The space motion of Ori OBI is mostly directed ra- 
dially away from the Sun, and the proper motion compo- 
nent is relatively small. The Hipparcos data did not allow 
de Zeeuw et al. ( |1999D to discriminate between the differ- 
ent subgroups in their selection procedure; they only give 
one proper motion and radial velocity for the whole Orion 
complex. It is possible that subgroup c has a motion that 
differs slightly from that of the other two subgroups, so 
that it cannot be ruled out as a candidate parent group. 
The age of subgroup c, ^5 Myr, is similar to the kinematic 
age of 53 Ari. By the argument given above this suggests 
that if Ori OBI c is the parent association of 53 Ari, then 
the formation mechanism is most likely the DES. 

In order to decide which of the Ori OBI subgroups 
is the parent group of 53 Ari, we need to know the dis- 
tances and velocities of the subgroups and the runaway 
star with a better accuracy than is now available. Fig- 
ure |l^ could then be used to pin down the parent group, 
and the mechanism which is responsible for the runaway 
nature of 53 Ari. Since subgroup a is the only one for 
which the BSS is indicated, finding a pulsar originating 
from subgroup a at the same time as 53 Ari would also 
clinch the issue. 
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Fig. 13. Predicted distance [left) and kinematic age 
(right) of ^ Per as a function of the distance of Per 0B2. 
The Hipparcos distance and its la error are indicated 
in the left panel. Distance estimates for Per 0B2 by 
Klochkova & Kopylov (1985, KK85), Cernis (1993, C93), 
and de Zeeuw et al. (1999, Z99) are indicated in the middle 
panel. 



Data: We adopt Wrad = 58.8 km s ^ for ^ Per (Bohan- 



nan & Garmany 1978; Garmany, Conti & Massey 1980 



Stone |1982| ; Gies & Bolton y^86|). This value differs by 
10 km s~^ from those quoted in the Hipparcos Input Cat- 
alogue (67.1 km s~\ Turon et al. [19 9 2D and the WEB 
catalogue (70.1 km s^^, Duflot et al. 1995), which derive 
from the value listed in the General Catalogue of Radial 
Velocities (70.1 km s"\ Wilson p^53| ) . We take the radial- 
velocity error to be 5 km s^^; this is equal to the ampli- 
tude of the velocity variations ind uced b y the non-radial 
pulsations of ^ Per (de Jong et al. 1999| ). The rotational 
velocity and helium abundance are Vi-ot sin i — 200 km s^^ 
and e = 0.18, respectively (see also Table S). 



5.2. ^ Persei & Perseus 0B2 



(star 3) 



The 07.5III star ^ Persei (HIP 18614) hes within the 
boundary of the Per 0B2 association on the sky. This 
positional coincidence, and the low density of early-type 
stars near Per 0B2, led Blaauw (1944) to propose ^ Per as 
a member of the association. At that time it was thought 
that the large radial- velocity difference between ^ Per and 
Per 0B2 (^--^40 km s~^) was due to uncertainties in the 
measurement of Uiad for ^ Per. However, when the radial 
velocity of ^ Per was confirmed, its membership of the 



Per 0B2 association became doubtful (Blaauw 1952a). In 
Paper I, Blaauw classified ^ Per as a runaway star from the 
Per 0B2 association, which naturally explains the discrep- 
ancy of the radial velocity of the star and the association. 
^ Per was the first star to be classified as runaway based 
on its Wrad alone. Most other runaways were recognized 
because their proper motions were directed away from an 
association. The parent group. Per 0B2 , has an age of 
^1 Myr ( e.g., Seyfert, Hardie & Grenchik 1960 ; de Zeeuw 
& Brand |1985| ). 



Simulations: Our orbit calculations (§|2) show that the 
kinematic age of ^ Per is ^1 Myr (Figure 13|). At that time 
the star was located ^5 pc from the center of Per 0B2, 
well inside the association. Figure |l^ also shows that the 
present distance of the runaway is 360 pc, assu ming 318 pc 
as the distance of Per 0B2 (de Zeeuw et al. |1999| ). This 
distance for ^ Per is consistent with the Hipparcos parallax 
at the 2(7 level. 

We infer that the BSS is responsible for the runaway 
nature of ^ Per based on (i) the 6 Myr age of Per 0B2 
at the time that ^ Per was ejected, (ii) the high helium 
abundance of ^ Per, (iii) its blue straggler nature (§^, 
and (iv) the large rotational velocity (see §^ . Further ev- 
idence of a supernova explosion in the Per 0B2 associa- 
tion is provided by a shell structure containing HI, dust. 



OH, CH, and other molecules (Sancisi 1970; Sancisi et al 



1974). This feature has been interpreted as a supernova 
shell which is physically connected to the Per 0B2 asso- 
ciation. We have not found a pulsar counterpart. 

^ Per presently illuminates the California Nebula 
(NGC 1499), resulting in an HII emission region. The dis- 
tance of this nebula is hard to determine (350-525 pc; 
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Fig. 14. Predicted distance iltfi) and kinematic age 
{right) of C Pup as a function of the distance of Tr 10. 
The Hipparcos distance and its la error are indicated in 
the left panel. The distance estimate by de Zeeuw et al. 
(1999, Z99) is indicated in the middle panel. 
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Fig. 15. Predicted distance {left) and kinematic age 
{right) of A Cep as a function of the distance of Cep 0B3. 
The Hipparcos distance and its la error are indicated 
in the left panel. The distance estimate by Crawford & 
Barnes (1970, CB70) is indicated in the middle panel. 



Bohnenste ngel fc Wendker 1976 ; Sargent 



1979; Klochkova 
1985|) , but must 



& Kopylov |1985| ; ShuU & van Steenberg 
be similar to that of ^ Per, i.e., ~360 pc. 



5.3. C Puppis, Vela R2, Vela 0B2 & Trumpler 10 

(star 10) 

The 041 star C Puppis (HIP 39429) is the brightest 
and nearest single O star to the Sun. Its location out- 
side any known association and its large space velocity 



(~60 km s ) led Upton (1971) to propose C Pup as 



a runaway star. In spite of the many groups of young, 
massive stars in the direction of C, Pup, Upton could not 
make a unique identification of the parent association. He 
proposed Vel 0B2 as the most likely candidate. In a re- 
cent s tudy of the stars and ISM in the direction of Vela, 
Sahu (1992) proposed the R association Vel R2 as a pos- 
sible parent for C, Pup, based on the retraced path of the 
runaway on the sky. For our simulations we considered 
Vel 0B2, Vel R2, the open cluster NGC 2 391, a nd the 
cluster/association Tr 10 (de Zeeuw et al. 1999) which 
also lies in the same direction as Vel R2. 

Data: Table || summarizes the data for C, Pup. The posi- 
tion and velocity for the associati ons V el 0B2 and Tr 10 
are adopted from de Zeeuw et al. ( 199S| ). Only two of the 
Vel R2 members are contained in the Hipparcos Cata- 
logue. Their distances, 411+^^| (HIP 43792) and 294+^f 
(HIP 43955), have large uncertainties or do not agree 
well wit h th e canonical distance of Vel R2, ^--^870 pc 
(Herbst |l97|). We are therefore unable to obtain mean- 
ingful phase-space coordinates of Vel R2. 

Simulations: Although we are unable to run our simula- 
tions for Vel R2, we conclude that it is not a likely candi- 
date parent association for Q Pup. The distance difference 
between Vel R2 (~870 pc) and C Pup (~400 pc, its canon- 
ical distance) is too large. The relative radial velocity be- 



tween the star and association is ~40 km s ^ (frad.c Pup = 
-23.9 km s^i and Wrad.vd R2 ^ 20 km s-i)|. The differ- 
ences in distance, ^400 pc, and velocity, ~40 km s~^, be- 
tween C Pup and Vel R2 yield a kinematic age of '--^lO Myr. 
This is older tha n the expected life time of C, Pup (van 
Rensbergen et al. 1996 ). Hence Vel R2 is not likely to be 



the parent group. 

Our simulations also show that Vel 0B2 is not the 
parent association. The minimum separation between the 
association and the runaway star is never smaller than 
40 pc for reasonable association distances. Since the as- 
sociation radius is, at maximum, 30 pc, we conclude that 
C Pup has never been inside the boundaries of Vel 0B2. 
We similarly rule out NGC 2391 as parent group. 

The simulations for the Trumpler 10 group result in 
minimum separations of >10 pc. The inferred kinematic 
age is ~2 Myr (Figure |lj). Ten parsec is comparable to 
the radius of Tr 10, so we cannot unambiguously iden- 
tify or exclude it as the parent association. Furthermore, 
if C Pup was in or near Tr 10, then its current distance 
must be 250-350 pc (Figure |lj) , which is smaller than the 
canonical distance of 400 pc. 

The Vela region contains many young stellar clusters, 
and suffers from a fair amount of extinction (although 
C Pup itself is almost unreddened). It is therefore reason- 
able to assume that we have not yet identified the parent 
group of C Pup. Similar conclusions were obtained by Van- 
beveren et al. (1998), and Vanbeveren, De Loore & Van 
Rensbergen (1998). 



5.^. A Cephei, Cepheus 0B2 & Cepheus 0B3 (star 22) 

A Cep (HIP 109556) was first classified as a runaway star 
by Blaauw (Paper I) . He noted that its position on the sky, 
the direction of the proper motion, and the radial velocity 



® We assume that Vel R2 has a motion similar to all 
other groups in the direction of the constellation Vela, i.e., 
~20 km s~^, almost independent of distance. 
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were consistent with an origin in the association Cep 0B2 
(see, e.g., figure 22 in de Zeeuw et al. 199E). This luminous 
supergiant (06I(n)f) is located below Cep 0B2 (in Galac- 
tic coordinates) and is traveling away from the Galactic 
plane and Cep 0B2. Several other stellar groups are lo- 
cated near Cep 0B2, and of these, Cep 0B3 may also be 
a possible parent of A Cep. We consider both associations. 

Data: The data for A Cep are given in Table ^. The 
phase-space co ordin ates of Cep 0B2 are adopted from de 
Zeeuw et al. ( 1999| ). Unfortunately, the Hipparcos data 
did not allow these authors to obtain meaningful results 
for Cep 0B3 which is at a distance of ~730 pc (Craw- 
ford & Barnes 1970 ). To estimate the phase-space coor- 
dinates of Cep 0B3 we used the mean position, proper 
motion, and radial velocity for th e Cep 0B3 members of 
Blaauw, Hiltner & Johnson ( |1959D : (£, b) = (110?50, 2?91); 
{ni*,Hb) = (-2.72±0.28,-1.78±0.30) mas yr"!; w^ad = 
-17.4±3.0 km s^^. 



objects, so our discussion is relatively brief. The results 
are summarized in Table ^. 



Simulations: The orbit calculations show that Cep 0B2 
cannot be the parent group of A Cep. The simulations 
do result in small minimum separations between the two, 
but only for Dx cop ~ 250 pc and a kinematic age of 
~9 Myr. These values do not agree with the observed pho- 



tometric (860 pc, Gies 1987) and trigonometric distances 



(505^95'^ pc, Hipparcos) of A Cep, nor with the age of the 



association (5-7 Myr, de Zeeuw et al. 1999| ). Moreover, 
A Cep is a massive O supergiant, and its lifetime cannot 
be more than a few million years. 

When we run the simulations using Cep 0B3 as the 
parent group we also obtain minimum separations <10 pc. 
Figure |l5| shows that the expected distance of the runaway 
is now ~45Q pc and that the kinematic age is ~4.5 Myr. 
This is a little on the large side for the nominal lifetime 
of a 40 Mq star, but might not be impossible. Cep 0B3 
consists of two subgroups with ages of 5.5 (subgroup b) 
and 7.5 Myr (s ubgroup a) (e.g., Jordi, TruUols & Galadi- 
Enn'quez 1996| ). Considering the high helium abundance 
and large rotational velocity of A Cep, subgroup a is a 
likelier parent of the runaway than subgroup b, since the 
age difference between the subgroup and the runaway is 
3 Myr for a. For subgroup b this is only 1 Myr, leaving 
little time for binary evolution. We conclude that A Cep 
is likely to have become a runaway star as the result of a 
supernova explosion in a binary system in subgroup a of 
Cep 0B3 ~4.5 Myr ago. 



6. Identification of nevif parent groups 

The orbit retracing technique allows us to identify the 
(likely) parent group for thirteen 'new' single runaways, 
one new pair, two more pulsars, and Geminga from the 
samples defined in 93. Little is known about most of these 



6.1. HIP 3881 & Lacerta OBI 



(star 1) 



The orbit of the spectroscopic binary v Andromedae 
(B5V-1-F8V) intersected the Lac OBI association 8- 
10 Myr ago. If this is the parent group, then the kinematic 
age is comparable to the (uncertain) age of the associa- 
tion ('--^lO Myr). This, together with a normal rotational 
velocity for the primary (80 km s~^) and the binary na- 
ture (two main-sequence stars) of this runaway, suggests 
that it was formed by dynamical ejection. 

6.2. HIP 22061, HIP 29678, Geminga & the A On region 
(stars 4 and 7 and neutron star 9) 

The two stars form an analog of the runaway pair AE Aur 
and ^ Col (§||). They also move in opposite directions and 
with similar space motions: 86.5 km s~^ for HIP 22061 
(B2.5V) and 63.0 km s^^ for HIP 296780 (BIV) (velocities 
are with respect to their own standard of rest). Retracing 
the orbits of HIP 22061 and HIP 29678, we find that the 
two stars were close together ^1.1 Myr ago, suggesting 
that these two runaways were also formed by the DES. We 
did not find a possible third party (i.e., another runaway 
in the neighbourhood or a massive binary). 

Applying the principle of conservation of linear mo- 
mentum at the time of the encounter, as we did in §|[ 
we can predict the properties of the parent cluster. We 
can only use these two stars and not three as in §4.5 



We find that the parent cluster should be located around 
(£, b) - (196?5, -12?0). This coincides wit h the A Orionis 
star- forming region (e.g., Gomez & Lada 1998; Dolan & 
Mathieu 1999D , which contains at least three young stel- 
lar clusters (the A Ori cluster and the clusters associated 
with the dark clouds B30 and B35) and is surrounded by 
the A Orionis ring. Several authors have suggested that 
this expanding ring of molecular clouds is the result of a 
super nova explosion ~0.35 Myr ago (e.g., Cunha & Smith 
1996 ). The predicted cluster position does not coincide 
with one of the three star clusters. Furthermore, the pre- 
dicted radial velocity of the cluster, ~10 km s~^, differs 
significantly from that of the A Ori clusters, ~24 km s~^. 
However, these differences might be erased if a third body 
(either a single star or a binary) was involved (cf. §|1)0. 

The conclusion that the DES is the acting mecha- 
nism for these runaways is supported by (i) the youth 



^ HIP 29678 is a visual double star. Gies & Bolton ( |1986| ) 
did not detect any radial velocity variations. 

Dolan & Mathieu (L999) report that long term monitoring 



of A Ori's radial velocity indicates that this star is a binary 
with a systematic velocity («rad ~ 18 km s~^) which differs 
by 6-8 km s~^ from that of the clusters in the A Ori region. 
However, the simulations show that A Ori cannot be related to 
the event that created the runaways HIP 22061 and HIP 29678. 



R. Hoogerwerf et al.: On the origin of the O and B-type stars with high velocities II 



23 



of the clusters in the A Ori star-forming region, 2-6 Myr 
(Dolan & Mathieu 1999| ), (ii) the density of these clusters 
(Dolan & Mathieu) , and (iii) the small rotational velocity 
of HIP 29678, <25 km s^^ (Morse, Mathieu, Levine |199lD . 

It is worth mentioning that HIP 22061 and HIP 29678 
are not the only objects running away from the A Orio- 



6.6. HIP 46950 & IC 2602 



(star 12) 



Somewhere between 2 Myr and 10 Myr ago the B1.5IV 
star occupied the same space as the open cluster IC 2602. 
The only bit of information available for the identification 
of the runaway's origin is the difference between the kine- 
matic age and the age of the open cluster (^25 Myr). This 
would suggest that the BSS is the origin. 



nis region. Frisch ( |1993D and Smith, Cunha & Plez ( [1994D 
suggested that the neutron star Geminga also originated 
from this star-forming region (Figure ^; but see Big- 

nami & Caraveo |l996| ). Moreover, the age of Geminga g 7_ jjjp j^ggj^g ^ i^wer Centaurus Crux (star 13) 



(^^350 000 yr) agrees well with the time of the supernova 
explosion which created the A Orionis ring 



6.3. HIP 38455 & Collinder 135 



(star 8) 



The ellipsoidal variable star HIP 38455 (B2V) moves away 
from the open cluster Collinder 135 with a velocity of al- 
most 40 km s^^, mostly in the radial direction. The orbits 
of the runaway and cluster intersected each other about 
3 Myr ago; this is significantly smaller than the age of the 
open cluster: ^--^35 Myr. This difference, and the large rota- 
tional velocity (wmt sini = 212 km s~^), suggest that the 
runaway originated via the BSS. Ellipsoidal variable stars 
show brightness variation because of their non-spherical 
shapes. The deformation of the star is thought to be due 
to tidal interactions with a companion star (Beech 19851 ). 
The fact that HIP 38455 is double (cf. Hafner & Drechsel 
1986| ) suggests that the kick velocity of the compact object 
created in the supernova explosion was small and that the 
binary remained bound. This system might in the future 
become a high-mass X-ray binary (sQ). 



6.4. HIP 38518 & Vela OB 2 



(star 9) 



This B0.5Iab supergiant is located behind the Vel 0B2 as- 
sociation and is currently moving away from it. Retracing 
the orbit shows that ~6 Myr ago the star was located in- 
side the Vel 0B2 association. This kinematic age is smaller 
than the association age (10-15 Myr). Adding to this the 
large rotational velocity (f rot sin i = 220 km s^i) and the 
blue straggler character of HIP 38518 (see §^ for a fur- 
ther discussion), we infer that the BSS is the preferred 
scenario. 

6.5. HIP 42038 & Upper Centaurus Lupus or IC 2391 

(star 11) 

The path of the B3V runaway star HIP 42038 traversed 
the open cluster IC 2391 ~6 Myr ago and the Upper Cen- 
taurus Lupus association ^8 Myr ago. Little is known 
about HIP 42038, making it difficult to determine its ori- 
gin. The only available criterion is that the kinematic age 
is smaller than both the age of IC 2391 (45 Myr) and the 
age of Upper Centaurus Lupus (13 Myr) favoring the BSS. 



HIP 48943 is a runaway from the Lower Centaurus Crux 
association. The orbit of the B5Ve star intersected the as- 
sociation center ^--^4 Myr ago. The runaway has a rotational 
velocity of 230 km s"^ and its kinematic age is significantly 
smaller than the age of Lower Centaurus Crux, ~10 Myr. 
The BSS is thus the most likely formation mechanism. 



6.8. HIP 49934 & IC 2391 or IC 2602 



(star 14) 



The retraced orbit of the emission-line star HIP 49934 
(B2IVnpe) intersects two open clusters: IC 2391 ^--^3 Myr 
ago and IC 2602 ~6 Myr ago. The large difference between 
the nuclear ages of both clusters (45 Myr and 25 Myr, 
respectively) and the kinematic age indicates that the BSS 
is the most likely origin of HIP 49934. This assumption 
is supported by the large rotational velocity (urot sin i = 
280 km s^i). 



6.9. HIP 57669 & IC 2602 



(star 15) 



About 3 Myr ago the emission-line star HIP 57669 left the 
open cluster IC 2602. This B3Ve star has a large rotational 
velocity (t^rotsini = 251 km s~^) and its kinematic age 
differs significantly from the age of IC 2602. Both these 
points suggest the BSS as the origin of HIP 57669. 

6.10. HIP 69491 & Upper Centaurus Lupus or 

Cepheus 0B6 (star 16) 

This candidate runaway is an eclipsing binary (B5V). Its 
path traverses the Upper Centaurus Lupus association and 
the Cep 0B6 association. The respective kinematic ages 
are 3 and 10 Myr. Both these kinematic ages are smaller 
than the association ages of 13 Myr for Upper Centau- 
rus Lupus and ~50 Myr for Cep 0B6. This discrepancy 
in ages combined with the binary nature of the candidate 
runaway excludes both the standard BSS and the DES 
and might suggest that neither Upper Centaurus Lupus 
nor Cep 0B6 is the parent of HIP 69491. The BSS is 
highly unlikely because the eclipses imply two objects of 
similar size and not, for example, a main-sequence star 
and a compact object. The DES is excluded because both 
associations were no longer compact/dense enough for dy- 
namical encounters to be efficient at the time of ejection. 
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What other mechanisms do exist to create a fast mov- 
ing (wspacc = 77 km s^^) binary system? One possibiUty 
is a supernova explosion in a triple system consisting of a 
hard binary and a third star with a larger semi-major axis 
(i.e., a stable triple system). This would result in either (i) 
a hard binary moving at moderate speed (<30 km s^^) or 
(ii) a fast runaway and a normal field star. In the latter 
case one of the stars in the binary explodes and creates 
a fast runaway. The third star, being weakly bound to 
the system would hardly be affected by the explosion. In 
the former case the single star explodes causing the bi- 
nary to start moving at the orbital speed it had within 
the triple system. This velocity should be small since the 
binary is much more massive than the third star. How- 
ever, neither case would create a runaway binary-system 
like HIP 69491. Whereas it is likely that the star origi- 
nated in either Upper Centaurus Lupus or Cep 0B6, the 
mechanism that formed this runaway remains unknown. 



6.11. HIP 76013 & Lower Centaurus Crux (star 17) 

Apodis is a Blnpe emission-line star, and is the bright- 
est component of a visual double system. This star has 
a radial velocity of 62 km and is moving away from 
the Galactic plane. Its orbit intersects the Lower Centau- 
rus Crux subgroup of Sco 0B2 2-3 Myr ago. Because of 
the large difference in kinematic age and association age 
(~10 Myr) the BSS is the most likely explanation for the 
runaway nature of Aps. 



6.14. HIP 102274 & Cepheus 0B2 



(star 21) 



6.12. HIP 82868 & IC 2602 



(star 19) 



Little is known about the BSVnpe star HIP 82868 whose 
orbit intersects that of the IC 2602 open cluster some 
6 Myr ago. A firm identification of the formation mech- 
anism is difficult, since we only know that the kinematic 
age differs from the age of IC 2602 (25 Myr). This suggests 
that HIP 82686 is a BSS runaway. 

6.13. HIP 91599 & Perseus 0B2 or Perseus 0B3 

(star 20) 

HIP 91599 is a known runaway star (B0.5V; Vitrichenko, 
Gershberg & Metik 1965); however, its parent associa- 
tion/cluster has never been identified. The simulations 
show that HIP 91599 originates from Per 0B2 or Per 0B3. 
The two associations have very different ages (^--^7 Myr for 
Per 0B2 and -50 Myr for Per 0B3). The kinematic age 
of HIP 91599 is ~8 Myr and -6 Myr for Per 0B2 and 
Per 0B3, respectively. Since we lack information on the 
rotational velocity and the helium abundance, we are un- 
able to conclude whether HIP 91599 is a DES runaway 
from Per 0B2 or a BSS runaway from Per 0B3. 



The B5 star HIP 102274 was located at the center of the 
Cep 0B2 association between two and three Myr ago, 
when the association was 3-4 Myr old. This kinematic 
age coincides with the time of the supernova explosion 
proposed by Kun, Balazs & Toth (1987) to explain the 
characteristics of the Cepheus bubble, a ring-like struc- 
ture of infrared emission. Taken together, this is strong 
circumstantial evidence for HIP 102274 being a BSS run- 
away. 

6.15. PSR J 0826+ 2637, PSR Jl 115+ 5030 & 

Perseus 0B3 (pulsars 1 and 4) 

The orbits of the pulsars J0826-h2637 and J1115-t-5030 
intersect that of the Per 0B3 association. This group con- 
tains the a Persei open cluster, and is often referred to as 
the Cassiopeia-Taurus association (de Zeeuw et al. 1999). 
The simulations show that if Per 0B3 is the parent of 
PSR J0826-I-2637, its kinematic age is ~1 Myr and its 
radial velocity is -100 km 8"^ For PSR J1115-f5030 we 
predict a kinematic age of ^--^1.5 Myr and a radial veloc- 
ity of —150 km s^^. The characteristic ages (P/(2P)) of 
the pulsars are 4.9 Myr and 10 Myr for PSR 10826-^2637 
and PSR J1115-I-5030, respectively. The unknown radial 
velocity of these pulsars makes it difficult to prove beyond 
doubt that these pulsars were born in the Per 0B3 associ- 
ation. Although Figure || shows that the orbits projected 
on the sky do not differ much for different radial veloci- 
ties, and that they cross Per 0B3, they may also cross the 
paths of other, more distant, associations or clusters not 
shown in the Figure. 

If the two pulsars orginated in Per 0B3, then the ini- 
tial periods would be 0.47 s for J0826-I-2637 and 1.53 s 
for J1115+5030, assuming no glitches occurred. The lat- 
ter value is large, which might indicate that this pulsar 
traveled longer, from another site of origin. 

It is not unlikely to find many pulsars associated with 
Per 0B3 since its age, —50 Myr, is comparable to the 
main-sequence life-time of an 8 Mq star. These are the 
least massive stars to explode as a supernova. Since the 
moment at which a star explodes, tsn, depends on its 
mass, (rsN ~ where a > and M > 8 Mq) and 

the number of stars of mass M, N{AI), also depends on 
the mass {N{M) - M"/^, where f3 > 0), the number of 
supernovae increases with time (A"sm ~ r'^/", for M > 
8 Mq). The number of supernovae, and thus the number 
of pulsars, will thus increase with time until the stars of 
8 Mq have exploded as supernovae. Afterwards the pulsar 
production rate will drop to almost zero. 



6.16. The Vela pulsar 



(pulsar 2 ) 



PSR J0835-4510 is only 10 000 yr old, and therefore has 
not traveled far from its birth place (—9' on the sky), the 
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Vela star-forming region at ~450 pc. It lies within the 
boundaries of the ~10 Myr old Vel 0B2 association (de 



8. Helium abundance versus rotational velocity 



Zeeuw et al. 1999), which is the likely parent group. 



7. Runaways and pulsars without parents? 

The sample of runaways and pul sars we have analysed 
here is severely incomplete (see §2J^). Furthermore, we 



have only a limited knowledge of where massive stars form 
in the Solar neighbourhood. The case of ( Pup described 
in §^ shows that some local parent groups may not yet 
have been found. Beyond ^500 pc, only the large com- 
plexes such as Cygnus 0B1/0B2 or Sco OBI are well 
documented. The accuracy of the three-dimensional posi- 
tions and velocities of these distant star-forming regions 
is poor, making determination of reliable orbits very dif- 
ficult. Runaway stars of spectral types BO and later are 
especially difficult to link to a parent group. These stars 
have main-sequence lifetimes of up to several tens of Myr 
and can therefore travel far (several kpc) from their places 
of origin. A similar argument holds for the pulsars. We 
suspect that most of the runaway stars and pulsars in our 
sample for which we were not able to identify a parent 
group originated outside the Solar neigbourhood. 

An example is the B2.5V star 72 Col, HIP 28756 (van 
Albada [l96]] the asterisk at [£,b] - [238°, -23°] in Fig- 
ure H). It has a peculiar velocity of ^200 km s^^ and its 
parent association is Sco OBI, at a distance of ^2 kpc 



(Humphreys 1978). Van Albada derived a kinematic age 
of ~14 Myr for 72 Col, based on a simple model of Galac- 



tic rotation (Kwee, MuUer & Westerhout 1954). This star 
does not appear in §§^ and ^ because its path did not 
carry it through one of the nearby associations. The Solar 
neighbourhood thus not only contains runaways for which 
the parent associations are also nearby, but it also contains 
runaways which originated far from the Sun. 

Another star that immediately catches the eye in Fig- 
ure | is HIP 94899 (the asterisk at , 6] - [341°, -26°]). 
This double star of spectral type B3Vn has a radial veloc- 
ity of 151 km s~^, and its path seems to cross the Per 0B3 
association. However, our simulations show that the run- 
away never comes within 40 pc of the association, implying 
that this system must have another, unknown, parent. 

We have seen in §| that PSR J1239-h2453 most likely 
originates outside the Solar neighbourhood. We did not 
find a parent group for the remaining pulsars because (i) 
an unreasonably large radial velocity (>500 km s~^) is 
necessary for the paths of the pulsar and parent group 
to intersect (PSR J1135-I-1551), just as we found for 
PSR J1239-f2453, or (ii) the past orbit simply does 
not intersect any of the nearby young stellar groups 
(PSR J0953+0755 and PSR J1456-6845). 



Blaauw ( 1993[ ) pointed out that most of the reliable run- 
aways known at the time have high helium abundances 
(e <; 0.13) and large rotational velocities (urotsini ^ 
150 km s^^), whereas most of the non-runaways have 
normal helium abundances and rotational velocities. Only 
one classical runaway star, AE Aur, does not follow this 
trend, and this was already suspected to be the result of 
the DES rather than the BSS. Blaauw suggested that the 
(e, Urot sini)-diagram of the small sample of known run- 
away stars supported the conjecture that most runaways 
are produced by the BSS, since these characteristics are 
natural con seque nces of close binary evolution (§0; van 
den Heuvel [l985| ). 

As shown in Table |^, only five of the 23 runaways listed 
there have a measurement of e, and six do not have a 
measured rotational velocity. In order to pursue Blaauw's 
suggestion, we therefore constructed a sampl e of O stars 
with known rotational velocities (Pe nny 1996 ) and helium 
abun dances (Kudritzki & Hummer 1990; Herrero et al. 
1992). We also determined, based on Hipparcos astrome- 
try and Hipparcos Input Catalogue radial velocities, the 
space velocities of these stars with respect to their local 
standard of rest. The (e, v^ot sin i) diagram in the left panel 
of Figure |l^ shows that these O stars can roughly be di- 
vided into three groups: (i) those with small rotational 
velocities, Uiot sini <, 80 km s^^, and normal helium abun- 
dances, e ^ 0.09, (ii) those with moderate rotational ve- 
locities, 80 Si Wiot sini ^ 150 km s^^, and normal to high 
helium abundances, 0.09 e 0.2, and (iii) those with 
large rotational velocities, Wrotsini ^ 150 km s^^, and 
high helium abundances, e <^ 0.14. The symbols in the left 
panel of Figure ^ are chosen according to the magnitude 
of the space velocity. The stars represented by filled circles 
have space velocities Wspacc > 30 km s^^, and the open cir- 
cles have Wspacc < 30 km s~^; the other symbols indicate 
stars for which no Hipparcos data are available (asterisks) 
or for stars with insignificant Hipparcos data (starred). 

The left panel of Figure |l^ does not show an obvi- 
ous separation between runaways (filled circles) and non- 
runaways (open circles). However, based on the available 
data for each star (position, radial velocity, distance mod- 
ulus, cluster membership) it is possible to decide whether 
the star is a runaway. The result is shown in the right panel 
of Figure |l^; it turns out that many of the stars with in- 
significant Hipparcos data are runaways (including some 
classical runaways), and that some of the stars indicated as 
high- velocity stars in the left panel are members of a clus- 
ter which has a peculiar motion. The resulting (e, Wrot sin i) 
diagram now shows a clear separation between the run- 
away stars and the normal O stars. Except for AE Aur 
(§^) all runaways have high helium abundances, e ^ 0.12, 
and large rotational velocities, Wiotsini ^ 150 km s~^. 
Only one star in this area of the (e, frotsini) diagram is 
not indicated as a runaway: the double star HIP 113306 
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^rot ^^'^ ' [km s~'] v^^^ sin i [km s~'] 

Fig. 16. Helium abundance (e: # relative to Hydrogen) versus rotational velocity for O stars. Left: open circles: stars 
with \v\ < 30 km s^^, filled circles: stars with \v\ > 30 km s^^, stars: stars with n — a-^ < mas, and asterisks: 
stars without Hipparcos data. Right: open circles: non-runaway stars, filled circles: runaway stars, and stars: doubtful 
runaways. The right of each panel display the helium abundances of stars with an unknown rotational velocity. See 
text for details. 



(0V7n). This star is a member of Cep 0B3 (e = 0.17, 
Urot sin I ~ 359 km s^^). 

We thus confirm Blaauw's conclusion that massive run- 
aways predominantly have high helium abundances and 
large rotational velocities, suggesting that they are formed 
mainly by the binary-supernova scenario. However, this 
conclusion is based on a limited sample which is by no 
means statistically complete. A systematic survey of the 
radial velocities, rotational velocities and chemical abun- 
dances of the early-type stars in the Solar neighbourhood 
is highly desirable. 

9. Blue stragglers 

One of the characteristics of runaway stars produced ac- 
cording to the binary-supernova scenario is that these 
stars are expected to be blue stragglers (§0). The mass 
transfer in close binary systems from the primary to the 
secondary prior to the supernova explosion deposits a 
large amount of hydrogen onto the future runaway. This 
new supply of fuel makes the runaway appear younger 
than the association/cluster in which it was born, i.e., it 
appears rejuvenated. 

Figure |l^ shows the colour vs. absolute magnitude di- 
agrams of the parent clusters discussed in this paper. The 
association members (dots) have been de-reddened follow- 
ing the Q-method; only the early-type members (AO and 
earlier) are shown. The solid lines denote the Schaller et 
al. ( |1992| ) isochrones for Solar metallicity and a standard 
mass loss rate for the ages of the associations. The run- 
away stars are denoted by starred symbols. A runaway can 
appear in more than one panel if two or more possible par- 



ents have been identified. The B — V colour and absolute 
magnitude have been determined using the spectral type 
of the runaways (Table |; Schmidt-Kaler |1982[ ). Three 
stars in Figure ^ clearly are blue stragglers: HIP 38518, 
Per, and A Cep; and three others could be blue strag- 
glers depending on the correct identification of the parent: 
C Pup, HIP 49934, and HIP 91599. The latter three stars 
have uncertain parent identifications (see §^ and §|^) . The 
blue straggler nature of the former three stars confirms 
their identification as BSS runaway (see Table |). The 
star C Oph has also been claimed to be a blue straggler 
(e.g., Blaauw 1993). However, the Upper Scorpius panel 
shows that ( Oph is the bluest star of the group, but it 
l ies on the main sequence, as also found by de Geus et al. 
( 1989| ) on the basis of uvhyjS photometry. 

By contrast to the BSS runaways, those produced by 
the DES are expected to follow the main sequence of the 
parent group. These runaways most likely did not experi- 
ence a period of binary evolution in which mass transfer 
was important. The runaway stars which we identified se- 
curely as DES runaways (AE Aur, ^ Col, HIP 22061, and 
HIP 29678) indeed fall on the main sequence of their par- 
ents (the Trapezium and the A Ori cluster). 



10. Concluding remarks 

We have used the Hipparcos astrometry for a sample of 
nearby candidate OB runaway stars to locate their parent 
groups, and to identify their formation mechanisms. We 
retraced the orbits of these candidate runaways, and de- 
termined where and when they passed through a possible 
parent group. We find that both mechanisms proposed for 
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Fig. 17. Colour vs. absolute magnitude diagrams of the runaways (stars) and their parent association/cluster (small 
dots). The association/cluster members have been de- reddened using the Q-method. The colours and absolute magni- 
tudes of the runaways have been determined using their spectral types (Table ||; Schmidt-Kaler 1983). The isochrones 
are from Schaller et al. (1992) for Solar metallicity and standard mass loss. The ages of the associations are indicated 
in the top right of each panel (US: Upper Scorpius, UCL: Upper Centaurus Lupus, LCC: Lower Centaurus Crux). 



the production of runaway stars, the binary-supernova sce- 
nario (BSS) and the dynamical ejection scenario (DES), 
operate. Table || summarizes the results. Even though the 
number of runaways discussed in this paper is small, and 
the weight of evidence varies, we find that roughly 2/3 of 
the runaways is produced by the BSS and 1/3 by the DES. 
This agrees with the results of binary population synthe- 
sis calculations by Portegies Zwart ( 20001 ). At present it is 
not feasible to extend this study to other runaways. This 
is mainly due to large uncertainties in the velocities and 
distances of the runaways, and the limited knowledge of 
star-forming regions with distances >500 pc. 

Tracing the runaway orbits back in time provides, for 
the first time, direct evidence that both scenarios pro- 



duce single runaway stars (Hoogerwerf et al. 2000). The 



orbit calculations demonstrate that the runaway C, Oph 
and the progenitor of PSR J1932-I-1059 once formed a bi- 
nary system in the Upper Scorpius association, and that 



the neutron star acquired a kick velocity of ~350 km s^^ 
in the supernova explosion. The runaways AE Aur and 
Col, and the binary l Ori were involved in a dynami- 
cal interaction (a binary-binary collision) ~2.5 Myr ago, 
which took place in the Trapezium cluster. 

The current investigation is biased towards finding 
BSS runaways. This is mainly due to the fact that the 
accuracy of the available data, and our knowledge of the 
location and motions of star- forming regions, restrict the 
study to ~700 pc. The small volume implies that we are 
only able to identify runaway stars with small kinematic 
ages of 0-10 Myr (i.e., runaways which recently left their 
parent association). Runaways which were created at an 
earlier time have most likely traveled outside our sam- 
ple limits. Since associations and open clusters can create 
BSS runaways during ~50 Myr (approximately the life- 
time of a 8 Mq star) and DES runaways only in the inital 
stages when the group still has a high density, we expect 
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Table 5. Parent associations, kinematic ages (tq) and for- 
mation mechanisms of the runaways and neutron stars dis- 
cussed in this paper. For runaway stars which have more 
than one possible parent we give the results for the individ- 
ual parent associations/clusters. The last column indicates 
the number of the runaway /pulsar in Figure |2| Abbrevia- 
tions: US: Upper Scorpius; UCL: Upper Centaurus Lupus; 
LCC: Lower Centaurus Crux; BSS: binary supernova sce- 
nario; DES: dynamical ejection scenario. 



HIP 


Name 


Parent 


To 
[Myr] 


Origin 


Fig.y 

LJ 


3881 


V And 


Laccrta OBI b 


9.0 


DES 


1 


14514 


53 Ari 


Orion OBI a 


4.3 


BSS 


2 






Orion OBI b 


4.8 


DES 


2 






Orion OBI c 


5.0 


DES 


2 


18614 


{ Per 


Perseus OB2 


1.0 


BSS 


3 


22061 




A Ori SFR 


1.1 


DES 


4 


24575 


AE Am- 


Trapezium 


2.5 


DES 


5 


27204 


^ Col 


Trapezium 


2.5 


DES 


6 


29678 




A Ori SFR 


1.1 


DES 


7 


38455 




Collinder 135 


3.0 


BSS 


8 


38518 




Vela OB2 


6.0 


BSS 


9 


39429 


C Pup 


7 






10 


42038 




UCL 


8.0 


BSS 


11 






IC 2391 


6.0 


BSS 


11 


46950 




IC 2602 


2-10 


BSS 


12 


48943 




LCC 


4.0 


BSS 


13 


49934 




IC 2391 


3.0 


BSS 


14 






IC 2602 


6.0 


BSS 


14 


57669 




IC 2602 


3.0 


BSS 


15 


69491 




UCL(?) 


3.0 


7 


16 






Cepheus OB6(?) 


10.0 


7 


16 


76013 


Aps 


LCC 


2.5 


BSS 


17 


81377 


C Oph 


US 


1.0 


BSS 


18 


82868 




IC 2602 


6.0 


BSS 


19 


91599 




Perseus OB2 


8.0 


DES 


20 






Perseus OB3 


6.0 


BSS 


20 


102274 




Cepheus OB2 


2.5 


BSS 


21 


109556 


A Cep 


Cepheus OB3 


4.5 


BSS 


22 




J0826+2637 


Perseus OB3 


1.0 


7 


1 




J0835-4510 


Vela OB2 


0.01 


7 


2 




J1115+5030 


Perseus OB3 


1.5 


7 


4 




J1932-I-1059 


US 


1.0 


BSS 


8 




Geminga 


A Orionis 


0.35 


7 


9 



to find more BSS than DES runaways because there are 
relatively many more old parents than young parents in 
the Solar neighbourhood. This bias is somewhat weakened 
by the fact that most dynamical interactions produce two 
runaway stars while the binary-supernova mechanism pro- 
duces only one. 

The creation of runaway stars modifies the mass func- 
tion of the parent group at the high-mass end, where the 
total number of stars is small. For example, the encounter 
in Orion described in §^ removed four stars with a total 
mass of order 70 Mq from the Trapezium cluster, while 
only six stars more massive than 10 Mq remain. Deriva- 
tion of the initial mass function of young stellar groups 
from the present-day mass function without accounting 
for the associated runaway stars leads to erroneous results. 

Our Hipparcos-based study has identified 56 runaway 
stars within 700 pc from the Sun, and tripled the subset 
of these for which a parent group is known (from 6 to 21). 
As mentioned in less than a third of the 0~B5 stars 
in the Hipparcos Catalog have a measured radial velocity. 



Obtaining these is likely to result in another factor of three 
increase in the size of the sample, so that statistical studies 
become possible. 

The next major step in our understanding of the origin 
of runaway stars will come when large datasets of micro- 
arcsecond (/xas) astrometry and accurate radial velocities 
(1-2 km s~^) become available. Distances accurate to a 
few parsec will allow for a final confirmation or rejection 
of the genetic link between runaways and their parents 
(e.g., Figure |lj). These data will become available over the 
next two decades with the launches of several astrometric 
satellites (FAME, SIM, GAIA). These aim to obtain /xas 
astrometry for a large number of stars, from 10 000 stars 
with SIM to 1 billion stars with GAIA. Besides astrometry, 
accurate radial velocities are also required; unfortunately, 
there is no dedicated effort to obtain these for a large 
number of O and B stars. 

The BSS and DES can produce runaway stars with 
spectral types beyond B5 (e.g., Kroupa 2000b; Portegies 
Zwart 2000). These will be harder to find, as the velocity 
distribution of the later-type stars in the Galactic disk is 
broader than for the 0~B5 stars, and the fractional pro- 
duction of low-mass runaways is small. Identifying their 
parent groups is also harder, because these stars may have 
traveled for much longer times. However, /ias accuracy 
astrometry complemented with accurate radial velocities 
will undoubtedly reveal such objects, and will provide fur- 
ther constraints on the binary fraction and the binary 
mass-ratios in open clusters and associations. 

Figure shows that there are 19 additional pulsars 
within one kpc for which an accurate proper motion is 
not available. A systematic program to measure these 
might allow the detection of more examples of pairs such 
as C Oph and PSR J1932+1059. It would also improve 
the characterisation of the pulsar population as a whole. 
VLBI techniques hold the promise of achieving sub-mas 
astrometry (positions, proper motions, and parallaxes) in 
the near future. 
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Appendix A: Distribution of differences 

In this appendix we record the distribution of the differ- 
ence between two observables, each of which has a Gaus- 
sian error distribution. We define two observable quan- 
tities, xi and X2, which are both distributed normally, 
Gi(xi; /ii, CTi) and G2{x2', ^J-2,o^2), respectively, where ni 
and (Ti are the mean and standard deviation of Gi, and 
A*2, 0-2 for G2- We define A = |xi — £2! as the absolute 
difference between the observables, where xi and X2 are 
the measurements of and X2- We calculate the distribu- 
tion of the difference, F{A), in the one- (i^iD(A)dA), two- 
(F2D(A)dA), and three-dimensional (F3D(A)dA) case. In 
all cases we took /ii = and fi2 = fi and ai = a2 = <J- The 
former simplification is harmless since only the difference 
/^i — M2 is important. The results are: 



Fid (A) = 



(see also tittp: //pulsar .princeton. edu/) 



exp 



exp 



(A + /i)^ ' 

4cr2 



(A.l) 
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Fig. A.l. Left to right: The distribution of absolute differences between two observable quantities, A, in one (Fid), 
two (f2D)j and three (J^sd) dimensions, taking into account the measurement errors. Top: ii — 0, I, 2, and 3 are 
represented by the solid, dotted, dashed, and dot-dashed lines, a is indicated in the top right of each panel. Bottom: 
(7 = 1,2, 3, and 4 are represented by the solid, dotted, dashed, and dot-dashed lines, /i is indicated in the top right of 
each panel. 



i^2D(A) = T^exp 

2(7^ 



1 A2_ 



^0 ( ^ ) , (A.2) 



where In is a modified Bessel function of order zero, and 



^^3d(A) 



exp 



exp 



1 (A - m)^ 

2 2ct2 

1(A + m)!" 

2 2cr2 
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Applying L'Hospital's rule we calculate the limit of i^3D 
for /i — > 0: 



lim F3D(A)dA 



„ _ o exp 



A^ 

'4a2 



(A.4) 



Figure Al and eqs ( |A.2| ) and (A.S) show that for the two- 
and three-dimensional cases there is a zero probability of 
measuring the same value for xi and X2, i.e., A = 0. 

Appendix B: The neutron star RX J185635-3754 



After this paper was submitted, Walter (2000) suggested 
that perhaps the past trajectory of the isolated neutron 
star RX J185635— 3754 intersected that of C Oph. He used 
HST WFPC2 observations over a three-year baseline to 
obtain a very accurate proper motion and parallax, and 
approximated the past orbit by a straight line. He con- 
cluded that the neutron star came very near C Oph 1.15 
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Fig. B.l. Left: Distribution of minimum separations, 
-Cimin(T-o), between C Oph and RX J185635-3754. The 
sohd hne denotes the expected distribution of Dmin, see 
j |3.3| . Right: Distribution of the times tq at which the min- 
imum separation was reached. 



Myr ago in Upper Scorpius for an assumed radial velocity 
of —45 km s^^. 

We have used the simulation machinery described in 
§H to analyse this case in the same way as done there 
for C Oph and PSR J1932+1059. We have again run 
three million simulations, covering the range —50 ± 50 
km s~^ for the radial velocity of the neutron star. Figure 
Bl presents the resulting distribution of minimum separa- 
tions r>,nin and associated kinematic ages tq. This figure 
can be compared directly with Figure We have seen in 
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§1 that 30822 simulations put PSR J1932+1059 within 
10 pc of C Oph, and in 4214 of these the encounter took 
place in Upper Scorpius. By contrast, only 748 simula- 
tions put RX J185635— 3754 within 10 pc of ( Oph, about 
1.5 Myr ago. None of these encounters occur within 15 pc 
of Upper Scorpius. We conclude that it is unlikely that 
RX J185635— 3754 is the remnant of the supernova that 
gave C Oph its large space velocity. We suspect that this 
neutron star formed long ago somewhere else in the Galac- 
tic plane. 



